DI7- 


SHIPBOARD  AND  SATELLITE  OBSERVATIONS  OF 

UPPER  OCEAN  VELOCITY  AND  TRANSPORT 

VARIABILITY  IN  THE  GULF  STREAM 


by 


3,adZ 


David  Michael  Schubert,  Jr. 

B.S.  United  States  Naval  Academy 
(1977) 

Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 

Master  of  Science 

at  the 

MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 

and  the 

WOODS  HOLE  OCEANOGRAPHIC  INSTITUTION 

September  1990 

©David  M.  Schubert,  Jr.  1990 

The  author  hereby  grants  to  MIT  and  to  WHOI  permission  to  reproduce  and 
to  distribute  copies  of  this  thesis  document  in  whole  or  in  part. 


^  \  Dr.  Carl  Wunsch 

Chairman,  Joint  Committee  for  Physical  Oceanography 


T248071 


GJL 


Table  of  Contents 


1.  Introduction 5 

2.  Description  of  Data 9 

2.1  Acoustic  Doppler  Current  Profiler  9 

2.2  Expendable  Bathythermograph  14 

3.  Estimation  of  the  Absolute  Velocity  Field    26 

3.1  ADCP  and  Geostrophic  Profiles  26 

3.2  Direct  Combination  of  ADCP  and  Geostrophic   Velocities  28 

3.3  Transport  Variability  42 

4.  Velocity  by  Combined  Inversion  Techniques    45 

4.1  The  Inverse  Problem  45 

4.2  The  Singular  Value  Decomposition  50 

4.3  Resolution  and  Variance  52 

4.4  Inverse  Results  for  Cruise  OC216  53 

4.5  Inverse  Results  for  Cmise  OC205  64 

5.  Temporal  Variability  between  Bermuda  and  Cape  Cod    73 

5.1  Gulf  Stream  Variability  73 

5.2  Ageostrophic  effects  76 

6.  Comparison  of  Sea  Surface  Height  and  Velocity  to  Satellite  Data 92 

7.  Conclusions    100 

2 


Shipboard  and  Satellite  Observations  of  Upper  Ocean  Velocity  and 
Transport  Variability  in  the  Gulf  Stream 

by 

David  Michael  Schubert,  Jr. 

Submitted  in  partial  fulfillment  of  the  requirements  for  the  degree  of 

Master  of  Science  in  Physical  Oceanography 

at  the  Massachusetts  Institute  of  Technology  and  the 

Woods  Hole  Oceanographic  Institution 

June  30,  1990 

Abstract 


Acoustic  doppler  velocities  are  combined  with  velocity  profiles  generated  from  XBT 
measurements  to  produce  estimates  of  the  flow  field  between  Bermuda  and  the  eastern 
coast  of  the  United  States.  Repeated  shipboard  measurements  along  an  ascending 
GEOSAT  subtrack  between  Bermuda  and  Cape  Cod  allow  study  of  rapid  Gulf  Stream 
Variability  along  the  track,  and  comparison  of  sea  surface  and  velocity  measurements  with 
those  computed  from  the  GEOSAT  altimeter.  The  shipboard  data  were  taken  during  two 
separate  cruises  on  the  R/V  Oceanus  in  April  and  December,  1989.  Using  mass 
conservation  constraints  and  inverse  techniques,  the  transport  across  the  Cape  Cod- 
Bermuda  track  has  been  balanced  with  transport  across  additional  ship  tracks  between 
Bermuda  and  Cape  Hatteras,  and  between  Bermuda  and  Nova  Scotia.  The  shipboard 
results  show  evidence  of  a  rapid  barotropic  mode  which  caused  changes  in  transport  along 
the  Cape  Cod-Bermuda  track  on  the  order  of  8  Sverdrups  in  a  week  period.  Comparisons 
of  sea  surface  velocity  and  dynamic  height  determined  from  the  ship's  data  with 
measurements  made  from  the  GEOSAT  altimeter  showed  a  consistent  picture  of  the  Gulf 
Stream  location  and  were  also  consistent  in  showing  smaller  scale  variations  in  flow.  The 
dynamic  height  difference  across  the  Gulf  Stream  was  approximately  10%  higher  for  the 
GEOSAT  measurements  than  for  the  shipboard  measurements,  which  is  within  the 
expected  errors  of  the  analysis  techniques. 
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Chapter  1 


Introduction 


Many  studies  have  been  made  over  the  years  of  water  transport  in  the  Gulf 
Stream.  Studies  done  in  the  1950's  and  earlier  used  the  dynamic  method,  which  measures 
the  horizontal  pressure  gradient,  and  uses  the  geostrophy  to  calculate  velocities  assuming 
the  force  caused  by  horizontal  pressure  gradients  in  the  water  is  balanced  by  Coriolis 
force.  Classically,  water  velocities  were  calculated  relative  to  some  deep  level  of  no 
motion.  Worthington  (1976)  analyzed  a  large  group  of  hydrographic  sections  across  the 
Gulf  Stream  to  investigate  the  seasonal  variation  in  Gulf  Stream  transport  relative  to  a 
level  of  no  motion  at  2000  meters.  It  is  now  well  known  that  significant  velocities  exist 
at  2000  meters  and  deeper,  making  estimates  such  as  those  by  Worthington  questionable 
as  representative  of  the  total  mass  flux.  In  the  1960's,  deep  current  measurements  were 
made  using  neutrally  buoyant  floats  in  connection  with  CTD  sections,  which  showed 
currents  of  speeds  between  6  and  11  cm/s  at  depths  of  approximately  2500  meters 
(Warren  and  Volkman,  1968).  Recentiy,  Halkin  and  Rossby  (1985)  have  made  a  study 
of  Gulf  Stream  transport  using  the  Pegasus  free-falling  velocity  profiler,  which  measures 
absolute  horizontal  velocity  profiles.  The  Pegasus  technique  is  quite  labor-intensive, 
requiring  an  ocean  bottom  transponder  array  to  acoustically  track  the  profiler,  and 
requiring  a  dedicated  research  vessel  to  deploy  the  profiler.  An  alternate  method  using 
hydrographic  data  and  an  Acoustic  Doppler  Current  Profiler  (ADCP)  for  direct  shallow 
water  velocity  measurements  was  reported  by  Joyce,  Wunsch,  and  Pierce  (1986)  and 


Pierce  and  Joyce  (1988).  This  method  also  requires  the  use  of  a  dedicated  research 
vessel,  but  does  not  require  as  much  ship  time  nor  any  external  transducer  network.  The 
work  presented  in  this  thesis  uses  a  variation  of  the  AE>CP/CTD/inverse  method  used  by 
Joyce,  Wunsch,  and  Pierce,  using  expendable  bathythermographs  (XBT's)  and  acoustic 
doppler  measurements  from  a  moving  research  vessel.  This  method  limits  measurements 
to  the  upper  ocean,  but  allows  faster  transits  to  be  made,  since  the  ship  does  not  need  to 
stop  for  CTD  stations.  Data  were  taken  along  a  GEOSAT  Altimeter  track  in  order  to 
investigate  the  possibility  of  using  direct  sea  surface  height  measurements  made  by 
satellite  to  provide  a  sea  surface  velocity  reference  for  geostrophically  determined  velocity 
profiles.  This  method  could  be  used  by  ships  of  opportunity  in  the  future,  eliminating 
the  need  for  a  dedicated  research  vessel. 

The  measurements  used  in  this  report  were  made  during  two  separate  cruises  on 
the  R/V  Oceanus.  The  first  cruise,  OC205,  conducted  in  April,  1989  consisted  of  two 
triangular  tracks  with  legs  running  from  Bermuda  to  the  200  meter  isobath  off  Cape 
Hatteras,  Cape  Cod,  and  Nova  Scotia  (Figure  1.1,  top)  The  central  leg  of  this  track  was 
chosen  to  coincide  with  an  ascending  subtrack  of  the  GEOSAT  altimeter.  This  leg  was 
traversed  three  times  during  the  April  cruise,  at  time  intervals  of  about  one  week.  The 
second,  somewhat  shorter  cruise,  OC216,  was  conducted  in  December,  1989,  and  covered 
the  same  area  as  the  lower  triangle  in  the  April  cruise  (Figure  1.1,  bottom).  During  the 
cruises,  ADCP  measurements  were  made  continuously  and  XBT's  were  launched 
approximately  once  an  hour.  The  GEOSAT  altimeter  is  in  an  orbit  which  causes  it  to 
repeat  its  track  over  the  earth  every  17  days.  During  the  April  cruise,  the  GEOSAT  pass 
occurred  between  the  second  and  third  ship  passage  along  the  track.  During  the 
December  cruise,  the  ship's  track  was  planned  so  that  the  Gulf  Stream  crossing  on  the 
final  leg  of  the  cruise  would  coincide  in  time  with  a  satellite  pass.  Unfortunately,  due  to 
equipment  malfunctions  associated  with  the  satellite's  age,  no  altimeter  data  were  received 
for  the  pass  during  the  OC216  cruise. 

The  geostrophic  velocity  profiles  generated  from  the  XBT  data  and  the  shallow 
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Figure  1.1  Tracks  for  R/V  Oceanus  cruises  in  April  1989  (top)  and  December  1989  (bottom) 


acoustic  doppler  velocity  profiles  will  be  combined  to  produce  an  absolute  velocity  field 
for  the  upper  760  meters  along  the  ship's  track.  The  direct  use  of  acoustic  doppler 
velocities  to  provide  a  reference  for  the  geostrophic  profiles  will  be  examined  first. 
Errors  in  the  doppler  velocities  using  this  direct  method  can  lead  to  substantial  errors  in 
total  transport.  The  geometry  of  the  sections  allows  us  to  impose  mass  budgeting 
constraints  on  the  transport  and  use  linear  inverse  theory  to  yield  a  more  accurate  picture 
of  the  flow  field. 

The  repeated  ship  legs  along  the  altimeter  track  allow  us  to  examine  variability 
along  this  track.  The  variability  section  uses  the  ADCP  measurements  to  estimate  the 
ageostrophic  components  in  the  momentum  equation,  and  show  how  the  Gulf  Stream 
varies  with  time. 

Finally,  the  sea  surface  dynamic  heights  and  velocities  computed  from  shipboard 
measurements  are   compared  to  those  computed  from  the  altimeter  in  order  to  validate 
the  use  of  altimeter  data  as  a  velocity  reference  for  hydrographic  data. 


Chapter  2 


Description  of  Data 


The  shipboard  data  consists  of  Acoustic  Doppler  Current  Profiler  (ADCP) 
velocities,  Expendable  Bathythermograph  (XBT)  temperature  profiles,  and  navigation 
positions  which  were  recorded  for  each  XBT  launch.  The  Acoustic  Doppler  system 
produces  profiles  of  horizontal  velocity  from  the  surface  down  to  a  depth  of  200  to  400 
meters.  The  XBT  temperature  profiles  are  combined  with  estimated  values  for  salinity 
to  produce  density  profiles,  from  which  geostrophic  velocity  profiles  can  be  obtained. 

Two  types  of  satellite  data  will  be  presented.  Sea  surface  height  measured  by  the 
GEOSAT  altimeter  will  be  differentiated  to  compute  sea  surface  velocities.  Infrared 
pictures  of  sea  surface  temperature  obtained  from  the  AVHRR  satellite  will  also  be  used 
qualitatively  to  examine  the  large  scale  features  of  the  Gulf  Stream  during  the  times  of 
the  other  observations. 

2.1    Acoustic  Doppler  Current  Profiler 

The  Acoustic  Doppler  system  computes  horizontal  velocity  components  by 
measuring  the  doppler  frequency  shift  of  sound  which  is  reflected  by  small  particles  in 
the  water.  Short  153.6  khz  pulses  are  transmitted  in  four  narrow  acoustic  beams,  which 
are  oriented  fore,  aft,  port,  and  starboard  of  the  ship.  Backscattered  sound  from  plankton 


and  other  small  particles  in  the  water  is  received  by  the  transducer  with  a  frequency  shift 
proportional  to  the  relative  velocity  between  the  ship  and  the  scatterers.  In  order  to 
compute  the  velocity  of  the  water,  the  ship's  velocity  must  be  accurately  known  and 
added  to  the  signal. 

The  ADCP  equipment  used  aboard  OCEANUS  is  the  RD  Instruments  Model 
VM0150  Current  Profiler.  In  order  to  reduce  the  effects  of  ship's  pitch  and  roll  on  the 
measurements,  the  pings,  which  have  a  pulse  repetition  frequency  of  one  second,  are 
averaged  over  a  two  minute  interval.  The  advertised  maximum  depth  for  the  ADCP 
system  operating  in  this  mode  is  350  meters.  During  the  cruises,  the  maximum  depth  at 
which  a  good  signal  (>25%  of  pings  good)  was  received  varied  between  200  and  400 
meters,  depending  on  the  roughness  of  the  weather  and  the  concentration  of  scatterers  in 
the  water. 

Water  velocities  are  determined  by  adding  the  ship's  velocity  to  the  doppler 
velocity.  Since  the  velocity  of  the  ship  is  typically  much  larger  than  the  velocity  of  the 
water,  errors  in  ship  velocity  can  lead  to  large  errors  in  the  water  velocity.  The  ship's 
velocity  is  determined  by  dividing  the  change  in  ship's  position  by  the  elapsed  time.  The 
accuracy  of  this  velocity  depends  on  precise  navigation.  During  the  cruises,  the  ship's 
position  was  determined  using  both  the  global  positioning  system  (GPS)  and  the  Loran-C 
navigation  system. 

In  order  to  compare  ADCP  velocities  with  geostrophic  velocities  computed  from 
the  XBT  data,  the  ADCP  data  were  vector  averaged  between  XBT  positions.  A  typical 
navigation  error  of  .2  km  between  2  XBT  positions  taken  one  hour  apart  would  cause  a 
5  cm/s  error  in  the  doppler  velocity.  An  error  in  XBT  position  will  introduce  equal  and 
opposite  velocity  errors  in  the  doppler  velocity  vectors  on  either  side  of  the  XBT.  To 
reduce  these  errors,  the  absolute  velocities  were  filtered  using  a  triangular  filter. 
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vn    =    -(vn_j  +  2vn  +  vn+1)  ,      for  all  but  end  velocities. 

(2.1) 

vLl    =    j(2v*nd  +  vend±i)  >      for  end  velocity  vectors 

The  effect  of  this  filter  is  to  reduce  the  navigational  noise  by  1/3  for  the  end 
velocity  vectors  and  by  half  for  all  the  rest  of  the  vectors.  This  filter  also  smoothes  the 
actual  velocities,  reducing  the  peak  and  broadening  the  width  of  the  Gulf  Stream. 


Another  source  of  error  in  the  ADCP  velocities  arises  from  misalignment  of  the 
transducer  with  the  ship's  gyrocompass.  A  calibration  run  (Joyce,  1988)  was  performed 
during  each  cruise  in  an  effort  to  eliminate  this  problem,  but  misalignment  effects  caused 
a  significant  bias  in  the  velocities,  especially  during  the  April  cruise.  With  the  ship 
moving  forward  at  a  speed  U,  a  small  error,  50  in  the  transducer  angle  will  generate  an 
apparent  velocity  perpendicular  to  the  ship's  track  of  Usin56.  This  error  will  be 
systematic,  appearing  as  a  consistent  velocity  error  to  one  side  of  the  ship.  A  transducer 
misalignment  of  .5°  with  a  typical  ship's  velocity  of  5  m/s  would  generate  a  cross-track 
velocity  error  of  4.4  cm/s. 

The  doppler  vectors  for  cruise  OC205  are  shown  in  Figure  2.1  In  this  Figure, 
only  the  second  of  the  three  central  track  legs  is  shown.  The  magnitude  and  direction  of 
the  velocity  vectors  along  the  200  meter  isobath  lead  one  to  suspect  a  systematic  bias. 
The  vectors  for  leg  2,  which  was  traveled  from  north  to  south,  point  consistently  off  the 
shelf,  while  the  vectors  for  leg  5,  which  was  traveled  from  south  to  north,  show  velocities 
onto  the  shelf.  In  each  case,  there  is  a  bias  to  the  left  of  the  ship's  track.  This  bias  will 
be  quantified  and  its  effects  removed  in  the  next  chapter.  The  Gulf  Stream  is  clearly  seen 
in  the  doppler  vectors  of  the  three  tracks  which  run  from  Bermuda  to  the  200  meter 
isobath.  A  large  cold  core  ring  is  also  observed  south  of  the  Gulf  Stream  on  leg  6 
between  Bermuda  and  Nova  Scotia. 


11 


0C205  LEGS  2-6  SMOOTHED  DOPPLER  VECTORS 


43    N 


39     N 


35     N- 


31    N4 


77    H 


73    H 


69    H 


65    M 


1.0  m/$ec 


Figure  2.1  Doppler  vectors  for  cruise  OC205  at  100  meters  depth.  The  vectors  for  the  Bermuda  to  Cape  Cod 
track  are  the  leg  4  vectors. 
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Figure  2.2  Doppler  vectors  for  cruise  OC216  at  100  meters  depth.  Legs  1  and  2  (top),  legs  3  and  4  (bottom). 
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The  doppler  vectors  for  the  December  cruise  are  shown  in  Figure  2.2.  Along  the 
Cape  Cod-Bermuda  satellite  track,  the  Gulf  Stream  is  about  100  km  further  to  the  North 
than  during  the  April  cruise.  There  are  two  large  cold  core  rings  near  the  track,  one 
centered  near  37°  N  which  moved  across  the  track  from  west  to  east  during  the  time 
between  legs  one  and  four,  and  one  west  of  Bermuda  which  remained  at  the  eastern  edge 
of  legs  two  and  three  during  the  cruise. 

2.2   Expendable  Bathythermograph 

Temperature  profiles  measured  by  expendable  bathythermograph  (XBT)  are  used 
to  determine  geostrophic  velocity  profiles.  Sippican  T-7  (765  meter)  XBT's  were  used 
on  both  cruises  except  along  the  shelf,  where  shallower  T-10  (200  meter)  XBT's  were 
used.  In  order  to  generate  density  profiles,  salinity  was  estimated  using  the  method  of 
Joyce,  et  al.  (1988).  The  procedure  for  generating  salinity  is  based  on  temperature  and 
depth: 

Starting  at  the  bottom  of  the  temperature  profile,  (765  m)  use  the  T/S 
relation  for  NW  Atlantic  Central  Water  (Armi  and  Bray,  1982)  to  estimate 
a  salinity. 

When  depth  reaches  200  meters,  hold  salinity  constant  to  the  surface 
unless  a  temperature  inversion  is  encountered. 

If  a  temperature  inversion  (>  0.5  °C)  is  encountered,  use  temperature  and 
depth  and  interpolate  towards  a  value  of  8  °C  /  32.5  psu,  which  is 
characteristic  of  shelf  water. 

This  technique  will  be  more  accurate  below  200  meters  than  above,  and  will  be 
more  accurate  in  the  Sargasso  Sea,  where  the  T/S  relation  is  accurate,  than  it  will  be  on 
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the  slope.   Joyce  et  al.  estimate  that  dynamic  height  errors  introduced  with  this  method 
will  be  of  order  2%  of  the  signal. 

Navigational  errors  will  also  affect  the  geostrophic  velocities,  but  not  to  as  great 
an  extent  as  the  ADCP  velocities.  The  same  .2  km  error  in  distance  between  XBT 
positions  described  above  will  produce  only  a  1.3%  error  in  velocity  for  XBT's  which  are 
15  km  apart. 

The  velocity  profiles  computed  from  the  XBT  data  assume  geostrophy,  and  do  not 
take  into  account  centripetal  acceleration  due  to  Gulf  Stream  curvature  or  any  other 
ageostrophic  components  in  the  momentum  balance.  These  components  can  be  significant 
and  are  discussed  in  chapters  5  and  6. 

Sections  of  temperature  and  potential  density  are  shown  in  Figure  2.3  through 
Figure  2.11.  All  of  these  sections  as  well  as  subsequent  sections  of  velocity  are  shown 
with  the  shore  end  of  the  track  to  the  left,  and  the  Bermuda  end  of  the  track  to  the  right. 
The  major  features  shown  in  the  doppler  vectors  are  also  apparent  in  the  temperature  and 
potential  density  sections.  Cold  core  rings  are  seen  in  the  OC205  sections  at  XBT  27-32 
of  leg  1,  XBT  118-123  of  leg  3,  and  XBT  254-258  of  leg  6.  There  are  also  cold  core 
rings  on  all  four  legs  of  the  OC216  cruise.  The  shelf  water  was  considerably  colder  near 
the  surface  during  the  April  cruise  with  surface  temperatures  ranging  from  8-15°  C 
compared  to  16-18°  C  for  the  December  cruise.  These  low  surface  temperatures  cause 
apparent  density  inversions  in  the  shallow  shelf  water  due  to  the  way  salinity  is  estimated 
in  the  upper  200  meters.  In  leg  4  of  cruise  OC205,  for  example,  the  potential  density 
section  shows  a  density  inversion  at  XBT's  190-191,  where  a,  at  the  surface  is  shown  to 
be  >  26.3.  The  actual  surface  value  of  a,  computed  using  surface  temperature  and  salinity 
measured  from  bottle  samples  is  26.1.  On  leg  7  of  cruise  OC205,  there  is  a  shallow  mass 
of  cold  water  between  XBT's  320-325  with  sea  surface  temperature  ~  7°  C.  Sea  surface 
potential  densities  computed  using  the  bottle  salinities  indicate  that  the  7°  C  water  is 
actually  less  dense  (a,  =  25.8)  than  the  surrounding  10°  C  water  (a,  =  26.0)  due  to  lower 
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Figure  2.3.    Temperature  and  potential  density  cruise  OC205  section  1.    Numbers  at  top  are  XBT  numbers. 
Distance  is  from  200  meter  isobath. 
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Figure  2.4.  Temperature  and  potential  density  cruise  OC205  section  3.   Distance  is  from  200  meter  isobath. 
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Figure  2.5.  Temperature  and  potential  density  cruise  OC205  section  4.  Distance  is  from  the  200  meter  isobath. 
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Figure  2.6.  Temperature  and  potential  density  cruise  OC205  section  6.   Distance  is  from  200  meter  isobath. 
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Figure  2.7.   Temperature  and  potential  density  cruise  OC205  section  7.   Distance  is  from  200  meter  isobath. 
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Figure  2.8.  Temperature  and  potential  density  cruise  OC216  section  1.  Distance  is  from  the  200  meter  isobath. 
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Figure  2.9.   Temperature  and  potential  density  cruise  OC216  leg  2.  Distance  is  from  the  200  meter  isobath. 
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Figure  2.10.  Temperature  and  potential  density  for  cruise  OC216  leg  3.  Distance  is  from  the  200  meter  isobath. 
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Figure  2.11.  Temperature  and  potential  density  cruise  OC216  leg  4.  Distance  is  from  the  200  meter  isobath. 
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salinity.  These  inversions  are  not  observed  in  the  OC216  sections  because  the  near- 
surface  temperature  gradients  are  smaller  and  the  errors  induced  by  the  salinity  estimation 
program  are  not  as  large.  In  the  next  section,  ADCP  velocity  profiles  will  be  compared 
to  geostrophic  velocity  profiles  for  several  XBT  pairs. 
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Chapter  3 


Estimation  of  the  Absolute  Velocity  Field 


3.1   ADCP  and  Geostrophic  Profiles 


The  geostrophic  velocity  profiles  were  generated  relative  to  an  assumed  level  of 
no  motion  at  765  meters.  Following  the  convention  of  Halkin  and  Rossby,  this  will  be 
referred  to  as  the  baroclinic  velocity  field,  and  the  765  meter  velocities  needed  to  obtain 
the  full  flow  will  be  referred  to  as  the  reference  velocity  field.  The  velocity  at  765  meters 
should  be  equal  to  the  difference  between  the  the  profile  of  ADCP  velocity  normal  to  the 
ship's  track  and  the  geostrophic  velocity  profile  generated  from  the  XBT  data.  With  no 
noise,  the  geostrophic  profile  would  have  the  same  shape  as  the  ADCP  profile,  offset  by 
the  actual  velocity  at  765  meters.  Errors  discussed  in  chapter  2  account  for  the 
differences  in  the  velocity  profiles.  Due  to  the  salinity  estimation  technique  used  with 
the  XBT  profiles,  the  geostrophic  velocities  are  likely  to  be  most  accurate  below  200 
meters  depth,  and  least  accurate  at  the  surface.  The  ADCP  velocities,  on  the  other  hand, 
are  most  accurate  above  200  meters,  with  velocity  dropping  to  zero  between  200  and  400 
meters  as  the  signal  is  lost.  The  reference  velocities  were  computed  using  an  average  of 
the  difference  between  ADCP  and  geostrophic  velocities  at  100,  150,  and  200  meters 
depth. 
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(b)  North  of  the  Gulf  stream,  cruise  OC205. 
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(e)  Cold-core  ring,  cruise  OC216  leg  1. 
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Figure  3.1.   Comparison  of  geostrophic  and  ADCP  profiles  of  velocity  normal  to  the  track.   Dots  are  ADCP 
velocities,  solid  black  line  is  geostrophic  profile. 


Velocity  profiles  for  several  station  pairs  are  shown  in  Figure  3.1.  The  top  two 
profiles  are  north  of  the  Gulf  Stream.  The  left  profile  is  from  leg  1  of  cruise  OC216  and 
shows  good  general  agreement  in  the  shapes  of  the  ADCP  and  the  geostrophic  profiles. 
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The  right  profile  is  for  XBT's  190-191  of  cruise  OC205,  where  the  apparent  density 
inversion  discussed  in  th  previous  chapter  is  observed  in  the  XBT  data.  We  can  see  that 
the  shapes  of  the  curves  are  quite  different  in  the  top  100  meters,  due  to  errors  in  the 
geostrophic  profile  induced  by  the  salinity  estimation  technique.  Profile  c  (mid  left)  is 
at  the  north  wall  of  the  Gulf  Stream,  and  shows  very  excellent  agreement  between  ADCP 
and  geostrophic  velocity.  Profile  d  (mid  right)  is  near  the  center  of  the  Gulf  Stream,  and 
shows  a  significant  difference  between  ADCP  and  geostrophic  velocity.  From  200  meters 
up,  the  geostrophic  velocities  are  decreasing  toward  the  surface,  while  the  ADCP 
velocities  are  increasing.  This  may  be  due  to  the  assumption  of  constant  salinity  above 
200  meters  used  in  computing  the  geostrophic  profile,  or  due  to  significant  ageostrophic 
terms  in  the  momentum  balance.  The  bottom  left  profile  (e)  is  in  a  cold  core  ring  south 
of  the  Gulf  Stream,  and  the  bottom  right  profile  (f)  is  from  the  Sargasso  Sea.  In  general, 
the  velocity  profiles  appear  consistent,  with  the  largest  differences  in  shape  occurring  near 
the  surface,  where  the  salinity  errors  are  the  greatest,  and  the  ageostrophic  terms  in  the 
momentum  equations  would  be  expected  to  be  largest. 

3.2  Direct  Combination  of  ADCP  and  Geostrophic 
Velocities 

The  straightforward  way  to  compute  the  total  velocity  field  once  the  reference  and 
baroclinic  components  are  calculated  is  simply  to  add  them  together.  Since  the  ADCP 
velocities  were  filtered  to  reduce  navigational  noise,  the  geostrophic  velocities  were  also 
filtered  before  adding  the  two  components  together.  Mass  transports  (M)  across  each 
track  are  computed  by  multiplying  the  velocities  (u)  by  the  cross  sectional  area  of  the 
track  (a)  and  by  the  density  of  the  water  (p).) 

M    =    uap  (3»1) 
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The  sections  were  set  up  to  be  able  to  take  advantage  of  mass  conservation 
constraints.  The  transport  across  the  sections  along  the  200  meter  isobath  should  be 
negligible  compared  to  the  transports  across  the  other  sections.  Beardsley  and  Boicourt 
(1981)  estimate  that  the  transport  over  the  whole  shelf  to  be  0.2  Sv.  Therefore,  the  net 
mass  transports  across  all  the  sections  running  from  the  shelf  to  Bermuda  should  be  the 
same.  When  the  transports  were  calculated  for  the  OC205  sections,  it  was  found  that  the 
transports  for  legs  3  and  6  were  much  larger  than  those  for  legs  1,  4,  and  7  (Table  3.1). 
Table  3.1   Cruise  OC205  section  transports.  Transports  are  in  kg/s  x  10^  =  Sverdrups. 


Section 

Baroclinic 
transport 

Total 
Transport 

After  bias 
removed 

1 

27.287 

26.885 

58.803 

3 

25.538 

95.752 

61.954 

4 

30.365 

22.183 

54.012 

6 

22.269 

82.632 

44.484 

7 

28.912 

22.745 

54.562 

The  large  imbalances  are  due  to  a  systematic  bias  in  the  doppler  velocity  discussed 
in  chapter  2.  Legs  1,  4,  and  7  were  traversed  from  South  to  North,  while  legs  3  and  6 
were  traversed  from  North  to  South.  A  bias  to  the  left  side  of  the  ship  would  add  to  the 
transports  on  legs  3  and  6  and  subtract  from  the  transports  of  legs  1,  4,  and  7.  The  value 
of  the  bias  velocity  which  most  nearly  balances  all  the  legs  is  4.48  cm/s  to  the  left  of  the 
ship.  When  this  bias  is  subtracted  from  the  reference  velocity  field,  the  transports  are  as 
given  in  the  last  column  of  Table  3.1.  The  resulting  velocity  sections  for  the  direct 
combination  of  ADCP  and  geostrophic  velocities  after  removing  the  bias  are  presented 
in  Figure  3.2  and  Figure  3.3. 
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Figure  3.2.  Velocity  normal  to  the  ship's  track  for  cruise  OC205  sections  1,  4,  and  7  by  direct  combination  of 
ADCP  and  geostrophic  data.  Distance  is  from  200  meter  curve.  Cross-hatched  areas  indicate  negative  velocity. 
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Figure  3.3.  Velocity  normal  to  ship's  track  for  cruise  OC205  sections  3  and  6  by  direct  combination  of  ADCP 
and  geostrophic  data.   Distance  is  from  200  meter  curve.   Cross  hatched  areas  indicate  negative  velocity. 
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OC205  Leg  1  Transport  and  Dynamic  Heights 
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OC205  Leg  4  Transport  and  Dynamic  Heights 
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OC205  Leg  7  Transport  and  Dynamic  Heights 
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Figure  3.4.   Transport  and  dynamic  heights  OC205  legs  1,  4,  and  7. 
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OC205  Leg  3  Transport  and  Dynamic  Heights 
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Figure  3.5.   Transport  and  dynamic  heights  OC205  legs  3  and  6. 


100 


33 


With  the  bias  removed,  significant  differences  in  transport  remain.  These  errors 
are  due  either  to  additional  errors  in  our  measurements,  or  to  short  term  variability  in  the 
flow  field.  In  order  to  examine  the  transport  of  each  section  more  fully,  it  is  necessary 
to  evaluate  the  transport  as  a  function  of  horizontal  distance  and  as  a  function  of  depth. 
Figure  3.4  and  Figure  3.5  show  the  transport  as  a  function  of  horizontal  distance,  along 
with  the  dynamic  height  computed  from  the  ADCP  velocities  at  100,  150,  and  200  meters 
depth  for  OC205  legs  1,  3,  4,  6,  and  7.  The  shape  of  the  transport  curves  closely  follow 
the  dynamic  height  curves,  with  small  differences  due  to  changes  in  transport  with  depth. 


The  transport  and  dynamic  height  curves  for  the  three  legs  between  Bermuda  and 
Cape  Cod  are  quite  consistent.  In  each  case  we  see  a  rise  in  dynamic  height  across  the 
Gulf  Stream  of  approximately  one  dynamic  meter  resulting  in  an  increase  in  transport 
across  the  Gulf  Stream  of  approximately  57  Sverdrups.  There  is  a  cold  core  ring  on  the 
track  during  leg  one,  which  causes  the  transport  to  dip,  then  rise  again  between  37°  and 
35°  N.  A  region  of  positive  velocity  between  33°  and  34°  N  causes  an  additional  rise 
in  transport  and  dynamic  height,  which  results  in  the  leg  1  transport  being  about  5 
Sverdrups  higher  than  those  of  legs  4  and  7.  The  slope  of  the  dynamic  height  across  the 
Gulf  Stream  in  leg  seven  is  less  than  for  the  other  legs,  because  the  Gulf  Stream  axis  has 
rotated  so  that  it  is  no  longer  nearly  perpendicular  to  the  ship's  track,  but  the  total  height 
difference  and  transport  is  consistent  with  the  other  two  legs. 

Comparing  the  transports  for  the  Cape  Cod-Bermuda  legs  to  the  other  two  Gulf 
Stream  crossings,  we  find  that  the  transport  for  the  Cape  Hatteras-Bermuda  leg  (leg  3)  is 
higher,  and  the  transport  for  the  Nova  Scotia-Bermuda  leg  (leg  6)  is  lower  than  the 
average  transport.  The  difference  in  transport  between  the  two  legs  appears  to  be  due  to 
effects  occurring  near  Bermuda.  At  approximately  100  km  from  Bermuda,  both  legs  3 
and  6  have  transports  of  about  55  Sverdrups  (Longitude  66°  W  for  leg  3,  latitude  33.5° 
N  for  leg  6  in  Figure  3.5).  From  those  points  to  Bermuda,  the  dynamic  height  and 
transport  for  leg  3  rises,  while  the  dynamic  height  and  transport  for  leg  6  drops  sharply. 
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This  appears  to  be  the  result  of  temporal  variability  in  the  currents  near  Bermuda.  The 
difference  in  time  between  the  turn  from  leg  3  to  leg  4  and  the  turn  from  leg  6  to  leg  7 
was  162  hours,  approximately  one  week.  Comparing  the  Bermuda  ends  of  legs  3  and  4, 
we  see  a  dip  then  a  rise  in  transport  coming  towards  Bermuda.  Legs  6  and  7  both  show 
a  steady  decrease  in  transport  south  of  34°  N,  although  the  effect  is  much  more 
pronounced  in  leg  6. 

The  transports  for  the  December  cruise,  OC216,  are  presented  in  Table  3.2.  The 
calibration  performed  for  this  cruise  appears  to  have  been  more  accurate  than  the  April 
cruise,  the  bias  which  rninirnizes  the  difference  in  transport  between  all  four  legs  being 
0.87  cm/s  to  the  left.  Again,  removing  the  bias  does  not  balance  the  transports,  leaving 
a  difference  of  10.5  Sverdrups  between  the  highest  transport  (leg  3)  and  the  lowest 
transport  (leg  1).  It  appears  that  something  happened  between  leg  2  and  leg  3  which 
appreciably  increased  the  transport  of  both  of  the  latter  two  legs  compared  to  legs  1  and 
2.  The  velocity  sections  for  cruise  OC216  after  removing  the  bias  are  presented  in 
Figure  3.6  and  Figure  3.7. 

Table  3.2  Cruise  OC216  Section  Transports.  Transports  are  in  kg/s  x  10*  =»  Sverdrups. 


Section 

Baroclinic 
Transport 

Total 
Transport 

After  bias 
removed 

1 

32.60 

59.40 

53.20 

2 

34.00 

49.49 

55.96 

3 

27.63 

70.28 

63.74 

4 

31.93 

54.96 

61.14 
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Figure  3.6.  Velocity  normal  to  ship's  track  for  cruise  OC216  sections  1  and  4  by  direct  combination  of  ADCP 
and  geostrophic  data.  Distance  is  from  the  200  meter  isobath.   Cross-hatched  areas  indicate  negative  velocity. 
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Figure  3.7.  Velocity  normal  to  ship's  track  for  cruise  OC216  sections  2  and  3  by  direct  combination  of  ADCP 
and  geostrophic  data.   Distance  is  from  the  200  meter  isobath.  Cross-hatched  areas  indicate  negative  velocity. 
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OC216  Leg  1  Transport  and  Dynamic  Heights 
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OC216  Leg  4  Transport  and  Dynamic  Heights 
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Figure  3.8.   Transport  and  dynamic  heights  OC216  legs  1  and  4. 
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OC216  Leg  2  Transport  and  Dynamic  Heights 
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OC216  Leg  3  Transport  and  Dynamic  Heights 
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Figure  3.9.   Transport  and  dynamic  heights  OC216  legs  2  and  3. 
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Plots  of  transport  and  dynamic  height  vs.  position  for  cruise  OC216  are  shown  in 
Figure  3.8  and  Figure  3.9.  The  Gulf  Stream  is  somewhat  stronger  than  during  the  April 
cruise,  with  a  transport  of  65  Sverdrups  and  a  dynamic  height  difference  at  100  meters 
depth  of  approximately  1.1  meter.  The  transport  and  height  profiles  for  the  two  Cape 
Cod-  Bermuda  legs  (legs  1  and  4)  are  very  consistent,  the  major  difference  being  a  steady 
drop  in  transport  on  leg  1  at  the  southern  end  of  the  track,  which  causes  the  total 
transport  to  be  8  Sverdrups  less  than  leg  4.  The  transports  and  dynamic  heights  for  the 
two  Bermuda-Cape  Hatteras  legs  show  quite  a  difference  in  transport  due  to  the  large  cold 
core  ring  between  32.5°  N  and  33.5°  N.  On  leg  3,  the  height  and  transport  lost  in  the 
northern  half  of  the  ring  are  returned  by  the  southern  half  of  the  ring,  as  would  be 
expected.  During  leg  2,  only  about  half  of  the  transport  removed  by  the  northern  half  of 
the  ring  is  returned.  This  appears  to  be  an  edge  effect.  The  ring's  cross-section  is  much 
wider  during  leg  2  than  leg  3,  and  its  center  is  further  to  the  east.  It  appears  that  during 
leg  2,  the  ship's  track  did  not  extend  to  the  end  of  the  ring.  By  leg  3,  the  ring  had  moved 
so  that  the  whole  width  of  the  ring  was  within  the  track,  and  the  transport  was  balanced 
between  the  northern  and  southern  halves  of  the  ring.  The  apparent  change  in  the  ring's 
diameter  is  most  likely  due  to  a  southward  movement  of  the  ring  such  that  leg  2  cut 
nearly  through  the  center  of  the  ring,  while  leg  3  crossed  the  ring  nearer  its  northern  edge. 

It  is  also  necessary  to  examine  the  vertical  structure  of  the  transport.  This  is  done 
by  dividing  the  sections  into  layers  defined  by  surfaces  of  constant  potential  density,  as 
described  by  Joyce  et.  al.  (1986).  If  we  assume  that  the  flow  is  along  surfaces  of  constant 
potential  density,  then  the  transport  within  each  layer  should  be  in  balance  as  well  as  the 
total  transport.  The  data  for  both  cruises  were  divided  into  six  layers  chosen  in  an  attempt 
to  resolve  the  major  water  mass  features.  These  layers  are  the  same  as  those  contoured 
in  the  density  sections  shown  in  chapter  2.  The  layer  transports  using  the  reference 
velocities  with  bias  corrections  for  both  cruises  are  listed  in  Table  3.3. 

We  can  see  from  the  table  that  there  are  significant  differences  between  layer 
transports,  as  well  as  the  total  transport  for  the  different  sections.    Although  a  small 
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Table  3.3  List  of  isopycnals  and  layer  transports.  Transport  units:  109  kg/s  *  106  m3/s. 


Individual  layer  transports  for  cruise  OC205 


Layer 
number 

<*e 

Leg  1 

Leg  3 

Leg  4 

Leg  6 

Leg  7 

1 

<26.3 

25.83 

24.91 

17.69 

14.57 

20.84 

2 

26.3-26.6 

14.37 

15.74 

10.16 

5.80 

10.77 

3 

26.6-26.9 

8.58 

10.28 

10.88 

8.93 

7.36 

4 

26.9-27.2 

9.33 

8.59 

8.59 

10.21 

9.52 

5 

27.2-27.6 

3.27 

2.14 

6.39 

6.55 

4.55 

6 

>27.6 

-2.58 

0.29 

0.31 

-1.58 

0.15 

Total 

58.80 

61.59 

54.01 

44.48 

54.56 

Individual  layer  transports  for  cruise  OC216 


Layer 
number 

<7e 

Leg  1 

Leg  2 

Leg  3 

Leg  4 

1 

<26.3 

22.00 

22.83 

22.33 

26.23 

2 

26.3-26.6 

9.24 

11.26 

14.59 

12.64 

3 

26.6-26.9 

8.36 

9.29 

11.65 

8.88 

4 

26.9-27.2 

6.86 

6.95 

8.95 

7.56 

5 

27.2-27.6 

5.31 

5.13 

5.38 

4.58 

6 

>27.6 

-1.44 

0.50 

0.83 

1.25 

Total 

53.32 

55.96 

63.74 

61.14 

amount  of  cross-isopycnal  flow  is  expected,  the  vertical  velocities  required  to  balance  the 
layer  flows  between  legs  are  too  large  to  be  physically  realistic.  In  order  to  better 
balance  the  total  and  individual  layer  transports,  inverse  techniques  will  be  used. 


41 


33   Transport  Variability 

Prior  to  using  individual  layer  and  total  section  mass  constraints  to  bring  the 
sections  closer  into  balance,  it  is  important  to  examine  why  the  section  transports  are 
different  and  how  much  change  would  be  expected  in  the  transports  in  the  time  required 
to  complete  the  sections.  Watts  (1983)  states  that  the  velocity  in  the  upper  800  meters 
in  the  Gulf  Stream  is  more  variable  than  is  apparent  looking  only  at  the  density  field  and 
concludes  that  this  variability  is  due  to  a  changing  barotropic  component  of  the  flow. 
From  Table  3.1  (page  29)  we  see  that  the  baroclinic  portion  of  the  transport  during  cruise 
OC205  varies  from  22.3  to  30.4  Sverdrups.  The  barotropic  part  of  the  transport,  which 
is  a  function  of  the  reference  velocities,  ranges  from  22.3  Sverdrups  (leg  6)  to  36.4 
Sverdrups  (leg  3),  a  range  of  over  14  Sverdrups,  which  is  almost  twice  as  much  as  the 
variation  in  the  baroclinic  transport.  There  does  not  appear  to  be  any  direct  correlation 
between  the  size  of  the  baroclinic  and  barotropic  components.  Leg  5  has  the  smallest 
amount  of  both  baroclinic  and  barotropic  transport,  but  leg  3,  which  has  the  second 
smallest  baroclinic  transport,  has  the  greatest  barotropic  component. 

Looking  at  the  repeated  tracks  between  Bermuda  and  Cape  Cod  (legs  1,  4,  and  7), 
we  can  see  how  much  of  this  variability  is  due  to  changes  in  the  Gulf  Stream  itself,  and 
how  much  is  due  to  changes  in  the  flow  field  outside  the  Gulf  Stream.  The  "Gulf  Stream 
transport"  for  legs  1,  4,  and  7,  determined  by  measuring  the  region  of  rapidly  increasing 
transport  (Figure  3.4),  is  59.4  Sverdrups  for  leg  1,  56.8  Sverdrups  for  leg  4,  and  57.0 
Sverdrups  for  leg  7.  Choosing  the  limits  of  the  Gulf  Stream  is  somewhat  arbitrary, 
especially  for  leg  7,  which  does  not  have  a  sharp  leveling  off  of  the  transport  increase 
south  of  the  Gulf  Stream.  The  transport  for  Leg  1  is  approximately  2.5  Sverdrups  higher 
in  the  Gulf  Stream,  and  approximately  5.5  Sverdrups  higher  overall  than  legs  4  and  7. 
Separating  the  transports  for  these  legs  into  baroclinic  and  barotropic  components,  the 
baroclinic  part  of  the  transport  varies  by  less  than  3  Sverdrups,  (Table  3.1)  whereas  the 
barotropic  part  varies  by  nearly  8  Sverdrups  between  leg  1  (31.5  Sv.)  and  leg  4  (23.6 
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Sv.).  This  seems  to  indicate  the  presence  of  some  large  scale  rapid  changes  which  could 
be  due  to  Rossby  waves. 

For  cruise  OC216,  we  can  compare  the  Gulf  Stream  transport  for  each  of  the  two 
sets  of  repeated  tracks  (Figures  3.8  and  3.9).  The  Cape  Cod-Bermuda  legs  have  Gulf 
Stream  transports  of  64.7  Sverdrups  (leg  1),  and  67.6  Sverdrups  (leg  4).  The  difference 
of  approximately  3  Sverdrups  between  legs  one  week  apart  is  similar  to  the  change  seen 
during  cruise  OC205.  The  Cape  Hatteras  tracks  have  Gulf  Stream  transports  of  67.1 
Sverdrups  (leg  2),  and  67.0  Sverdrups  (leg  3).  The  transports  are  nearly  identical,  which 
is  not  surprising  considering  the  small  elapsed  time  between  the  two  tracks.  (The  two 
Gulf  Stream  crossings  were  completed  in  just  over  one  day.) 

The  major  cause  of  variability  in  the  OC216  legs  is  caused  by  the  cold  core  ring 
near  Bermuda,  which  which  moved  to  the  south-west  during  the  cruise.  The  average 
transport  for  legs  3  and  4  is  approximately  8  Sverdrups  higher  than  for  legs  1  and  2  due 
to  edge  effects  caused  by  the  movement  of  this  ring  between  legs  2  and  3.  The  change 
in  transport  due  to  this  edge  effect  is  barotropic.  In  fact,  the  baroclinic  part  of  the 
transport  is  somewhat  smaller  for  legs  3  and  4  than  it  is  for  legs  1  and  2  (Table  3.2,  page 
35). 

In  summary,  the  following  types  of  variability  are  seen  in  the  data: 

Variability  in  Gulf  Stream  transport  between  Bermuda  and  Cape  Cod  of  order  3 
Sverdrups  in  a  one  week  period. 

Variability  due  to  edge  effects  of  approximately  8  Sverdrups,  caused  by  the  motion 
of  a  cold  core  ring  during  cruise  OC216. 

Large  scale  barotropic  variability  between  Cape  Cod  and  Bermuda  of  order  8 
Sverdrups  in  one  week. 
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The  inverse  models  of  the  next  chapter  will  use  mass  balance  constraints  to 
balance  the  flow  between  sections.  The  sources  of  variability  listed  above  must  be  taken 
into  consideration  when  making  the  models.  For  cruise  OC216,  there  is  a  large  change 
in  transport  between  legs  2  and  3  due  to  edge  effects  which  should  not  be  removed  by 
attempting  to  balance  all  four  legs  together.  Therefore,  the  inverse  model  for  cruise 
OC216  will  involve  only  balancing  leg  1  with  leg  2  and  leg  3  with  leg  4.  If  we  examine 
the  legs  two  at  a  time,  the  transport  differences  between  legs  1  and  2  are  small,  as  are 
the  differences  between  legs  3  and  4.  We  can  use  a  fairly  tight  constraint  on  the  total 
transport  when  conducting  the  two  leg  balances. 

For  cruise  OC205,  a  three  leg  balance  will  be  conducted,  applying  mass  constraints 
to  balance  the  flow  in  legs  3,  4,  and  6.  The  time  from  the  start  of  leg  3  to  the  end  of  leg 
6  was  9.5  days,  while  the  average  time  for  legs  1-2  and  legs  3-4  of  cruise  OC216  was 
less  than  six  days.  Due  to  the  much  larger  initial  imbalances  in  transport  for  the  three 
legs,  as  well  as  the  longer  period  of  time  involved,  less  weight  will  be  given  to  the  total 
mass  balance  constraints  when  constructing  the  inverse  model  for  cruise  OC205. 
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Chapter  4 


Velocity  by  Combined  Inversion  Techniques 


4.1   The  Inverse  Problem 


The  inversion  technique  combines  the  acoustic  doppler  velocities  with  the 
geostrophic  profiles  while  at  the  same  time  attempting  to  satisfy  certain  other  physical 
constraints  within  some  estimated  errors.  All  the  constraint  equations  are  set  up  in  a 
linear  matrix,  which  is  solved  using  least  squares  techniques,  specifically  the  Singular 
Value  Decomposition  (SVD),  which  yields  information  on  the  resolution  of  each  equation 
and  unknown.  These  inversion  procedures  were  first  applied  to  oceanographic  data  by 
Wunsch  (1978).  The  procedures  here  closely  follow  Joyce,  et.  al.  (1986),  and  Pierce  and 
Joyce  (1988)  who  applied  the  inverse  technique  to  a  combination  of  ADCP  and 
hydrographic  data.  There  are  several  differences  between  the  current  model  and  previous 
models.  Due  to  the  limits  on  XBT  depth,  only  the  transport  in  the  upper  765  meters  is 
being  balanced.  Also,  the  bias  levels  discussed  in  chapter  3  have  been  removed  to  make 
an  initial  correction  to  the  total  transports  prior  to  performing  the  inverse  calculations. 

The  model  parameter  matrix  is  made  up  of  equations  chosen  to  satisfy  the 
following  physical  constraints: 
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The  760  meter  velocities  are  consistent  with  those  determined  by  the 
differences  between  the  ADCP  and  geostrophic  velocity  profiles  and 
corrected  for  bias  errors. 

The  total  mass  transport  in  the  upper  760  meters  is  conserved  for  the 
triangular  areas  between  Bermuda  and  the  shelf. 

Mass  is  conserved  within  each  density  layer. 

The  model  is  written  in  matrix  format  as 

Ax    =    b    ,  (4.1) 

where  A  is  the  matrix  of  model  parameters,  b  is  a  vector  of  observations,  and  x  is  the 
vector  containing  the  unknowns.  (Throughout  the  text,  vector  quantities  will  be  denoted 
by  boldface  small  letters,  and  matrices  will  be  denoted  by  boldface  capital  letters.)  The 
unknowns  in  our  model  will  be  the  horizontal  reference  (765  m)  velocities  for  each 
station  pair,  and  the  cross-isopycnal  velocities  associated  with  the  chosen  density  surfaces. 
The  number  of  unknowns,  therefore,  will  be  equal  to  the  number  of  XBT  pairs  in  the  legs 
being  balanced,  plus  the  number  of  isopycnal  interfaces.  The  concept  of  cross-isopycnal 
flow  is  discussed  in  Wunsch  et.  al.  (1983).  The  horizontal  reference  velocities  will  be 
denoted  as  Xj  where  j  =  1  to  nsta,  the  number  of  station  pairs  in  the  legs  being  balanced. 
The  cross-isopycnal  velocities  will  be  denoted  by  w*,  where  i  is  the  number  of  the 
interface  between  density  layers.  Since  we  have  divided  the  sections  into  six  density 
layers,  there  will  be  five  interfaces  for  each  area  used  in  the  transport  balance. 

The    condition   that   reference   level   velocity   is   consistent   with   the    ADCP 
measurements  is  expressed  as 
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a j    =    qj    ±    Ej  ,       ;*  =  1,  nsta 


(4.2) 


where  cjj  is  the  true  reference  velocity,  o^  is  the  initial  estimate  of  reference  velocity  as 
determined  in  the  previous  chapter,  and  ei  represents  the  error  in  the  acoustic  doppler 
velocity.  As  discussed  previously,  the  expected  random  error  in  the  acoustic  doppler 
velocities  is  ±5  cm/s. 

The  transport  constraints  are  expressed  as  equations  balancing  the  inflow  and  the 
outflow  of  mass  in  any  given  layer,  or  the  whole  section. 


nsta 


(4.3) 


E  ay%  +  V  +  hiwt    ~  hi-iwi-i  =  °    • 


We  define  a,j  as  the  vertical  area  of  layer  i  within  station  pair  j  multiplied  by  the  density 
of  layer  i,  and  hj  as  the  horizontal  area  of  interface  i  between  two  ship's  tracks  multiplied 
by  the  density.  The  quantity  (cy  +  vt^)  is  the  sum  of  the  reference  and  baroclinic  velocity 
components  in  area  ij.  There  will  be  one  mass  conservation  equation  for  each  layer,  and 
one  for  the  total  section  mass  transport.   We  can  express  this  set  of  equations  as 


Axx    -    bx    «    0 


(4.4) 


where  A1  is  a  matrix  made  up  of  the  elements  a,j  and  h,,  -bj  is  the  initial  imbalance 
between  layer  transports  due  to  the  baroclinic  component  of  velocity,  and  x  is  the  set  of 
unknowns 


(4.5) 


x    = 


w 
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The  approximately  equal  sign  in  equation  (4.4)  represents  the  fact  that  mass  is  conserved 
only  to  within  some  level  of  error.  Equation  (4.4)  represents  the  pure  geostrophic  inverse 
problem,  where  property  conservation  properties  are  used  to  solve  for  reference  velocities. 
The  set  of  equations  (4.2)  can  be  added  to  the  matrix,  combining  the  ADCP  constraints 
with  the  mass  balance  constraints.  The  problem  can  now  be  stated: 


',«.=» 


\ 


x    = 


a 


b, 


(4.6) 


or  simply  Ax  =  b.    Solving  for  the  unknowns  is  called  the  inverse  problem,  because  it 
involves  inverting  the  matrix  A. 


Prior  to  solving  the  inverse  problem,  the  equations  must  be  weighted.  A  row 
weighting  factor  is  introduced  for  each  equation  which  is  inversely  proportional  to  the 
estimated  error  in  the  constraint  for  that  equation.  When  all  the  equations  are  multiplied 
by  their  weighting  factors,  the  expected  error  in  each  equation  is  normalized.  The 
equations  are  now  non-dimensional.  Column  weighting  is  also  imposed  on  the  A  matrix 
in  the  interest  of  numerical  stability.  There  is  an  artificial  tendency  for  the  solution  to 
favor  the  terms  in  the  equation  with  the  largest  magnitudes.  Therefore,  any  station  pairs 
which  were  closely  spaced  or  shallow  would  be  less  resolved  than  those  spaced  farther 
apart  or  deeper.  Additionally,  the  horizontal  interface  areas  are  orders  of  magnitude 
larger  than  the  vertical  areas,  and  the  horizontal  velocities  q^  are  much  larger  than  the 
vertical  velocities  w*.   The  column  weighting  scales  the  size  of  each  unknown. 

Let  the  row  weighting  matrix  be  R,  a  diagonal  matrix  with  elements  r^  =  aj5  the 
expected  error  in  the  jth  equation.   With  row  weighting,  equation  (4.1)  becomes 
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RAx    =    Rb    .  (4-8) 


To  weight  the  columns,  define  another  diagonal  matrix  C,  such  that  the  diagonal  elements 
c„  are  equal  to  the  area  of  the  station  pair  or  layer  associated  with  column  i.  C  has  a 
square  root  C*,  and  C^C*  =  I,  the  identity  matrix.  With  column  and  row  weighting, 
equation  (4.1)  becomes 

--    -  (4.9) 

RAC  2C2x    =    Rb    ,  y     } 


which  we  can  rewrite  in  the  form  of  equation  (4.1)  as 

A'x'    =    b1   ,  <4-10> 

where  A;  =  RAC'*,  x;  =  C^x,  and  W  =  Rb.  In  the  subsequent  discussion  the  primes 
will  be  dropped,  and  the  unprimed  values  will  refer  to  the  weighted  matrix. 

Our  system  will  generally  have  more  equations  than  unknowns,  therefore  the 
problem  is  formally  overdetermined.  The  conventional  method  for  solving  an 
overdetermined  problem  is  to  use  the  least  squares  approach.  The  least  squares  solution 
is 

i    =    [ATAVlATb    ,  t4-11) 

where  AT  is  the  transpose  of  the  data  matrix  A.  The  least  squares  solution  gives  the  best 
estimate  of  the  solution  in  the  sense  that  the  sum  of  the  square  differences  between  Ax 
and  b  is  minimized.  Although  our  A  matrix  is  formally  overdetermined,  the  equations 
in  the  matrix  may  not  be  independent.    If  the  number  of  independent  equations  is  less 
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than  the  number  of  unknowns,  the  problem  is  underdetermined.  In  solving  the  system, 
we  use  a  solution  technique  which  is  equivalent  to  the  least  squares  technique,  but  is 
equally  adept  at  handling  the  overdetermined  and  underdetermined  cases. 

4.2   The  Singular  Value  Decomposition 

The  Singular  Value  Decomposition  (SVD)  of  a  matrix  is  a  very  powerful  tool  for 
evaluating  the  resolution  and  variance  of  the  solution  when  solving  the  inverse  problem. 
The  Singular  Value  Decomposition  decomposes  the  matrix  A  into  two  sets  of  orthonormal 
eigenvectors  and  a  set  of  eigenvalues.  The  value  of  this  decomposition  is  the  ease  with 
which  it  lends  itself  to  qualitative  ranking  of  the  information  content  of  the  system.  The 
SVD  solution  simultaneously  minimizes  the  misfit  in  each  equation  and  the  deviation 
from  the  initial  model.  In  this  case,  the  magnitude  of  the  reference  velocity  corrections 
and  of  the  cross-isopycnal  velocities  is  rninimized.  Thorough  discussions  of  the  Singular 
Value  Decomposition  are  found  in  Lanczos  (1961),  Wiggins  (1972),  and  Wunsch  (1978). 

The  Singular  Value  Decomposition  of  A  can  be  written  as 

(4.12) 


A 

U 

A 

VT 

m  xn 

m  xm 

m  xn 

n  xn 

where  U  is  composed  of  m  orthonormal  column  vectors  u,,  and  V  is  composed  of  n 
orthonormal  columnvectors  v,.  The  matrix  A  is  a  diagonal  matrix  with  elements  \,  called 
the  singular  values  of  the  matrix.  The  U  and  V  matrices  form  a  coupled  eigenvalue 
problem: 

AATU=k2U,     ATAV  =  \2V    .  <4-13> 
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In  general  the  number  of  rows  m  is  not  equal  to  the  number  of  columns  n.  The 
number  of  non-zero  singular  values,  k  is  known  as  the  rank  of  the  matrix  A  and 
represents  the  number  of  independent  equations  in  the  matrix.  The  singular  values,  which 
are  all  positive,  are  arranged  from  largest  to  smallest  such  that  \+l  <  \.  The  rank  of  the 
matrix  is  always  less  than  or  equal  to  the  smaller  of  m  and  n.  As  will  be  shown  below, 
the  smallest  singular  values  can  introduce  large  variance  into  the  solution.  This  effect  can 
be  eliminated  by  setting  to  zero  all  singular  values  smaller  than  some  cutoff,  thereby 
reducing  the  rank  of  the  system. 

The  SVD  solution  of  the  inverse  problem  is  given  by 

(4.14) 
£    =     VA1UTb    . 


This  solution  is  equivalent  to  the  least  squares  solution,  and  also  satisfies  the  constraint 
that  it  minimizes  the  sum  of  the  squares  of  the  solution  parameters  (i.e.,  xTx  is  as  small 
as  possible).  Equation  (4.14)  can  be  rewritten  as  a  sum  of  the  eigenvectors 


J^        Ul    'b  tA    Iff. 

t    =    £  _L_V(    .  (4.15) 


From  this  expression,  it  is  easy  to  see  that  the  smallest  values  of  A,  will  have  the  largest 
effects  on  the  solution.  Choosing  the  cutoff  below  which  the  eigenvalues  will  not  be  kept 
in  the  summation  is  the  problem  of  determining  the  rank  of  the  system.  Wiggins  (1972) 
and  Lawson  and  Hanson  (1974)  discuss  several  approaches  to  making  this  decision.  As 
more  and  more  eigenvalues  are  kept,  the  residual  norm  HAx  -  bll  decreases  while  the 
solution  norm  llxTxll  increases.  The  best  choice  of  rank  is  one  for  which  the  residuals 
have  been  made  acceptably  small,  but  the  solution  has  not  been  made  too  large. 
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4.3   Resolution  and  Variance 

The  U  and  V  matrices  provide  useful  information  about  the  contributions  from  the 
equations  and  the  resolution  of  the  unknowns.   Since  the  eigenvectors  are  orthogonal, 

UUT  =  /_    ,  VVT  =  /.    ,  (4.16) 


m      »  n 


where  Im  is  the  m  x  m  identity  matrix,  and  I0  is  the  n  x  n  identity  matrix.  If  the  rank  k  is 
less  than  m  or  n,  only  the  first  k  columns  of  U  and  V  are  used  by  the  SVD.  If  we  denote 
the  matrices  made  of  the  first  k  columns  of  U  and  V  as  Uk  and  Vk,  we  find  thatUfcUj^I^ 
and  VkVj*In.  It  can  be  shown  that  the  diagonal  elements  of  VkVj  give  the  resolution  of 
the  n  unknowns  of  the  model.  At  rank  ny  all  the  diagonal  elements  of  VkVk  are  1  and  all 
the  unknowns  are  perfectly  resolved.  At  a  lower  rank,  some  of  the  diagonal  elements  of 
VkVk  will  be  less  than  1  and  the  corresponding  unknowns  will  be  less  than  perfectly 
resolved.   The  sum  of  the  diagonal  elements  of  VkVk  =  k. 

The  UkUk  matrix  gives  similar  information  about  the  contributions  of  each  of  the 
equations  to  the  solution.  The  m  diagonal  elements  of  UkUk  give  the  contribution  of  the 
m  equations  to  the  solution.  A  diagonal  element  of  1  means  the  corresponding  equation 
contributes  fully  to  the  solution.  A  diagonal  element  of  UkUk  near  zero  means  that  the 
corresponding  equation  is  not  independent  of  the  other  equations  and  is  not  contributing 
to  the  solution.   The  sum  of  the  diagonal  elements  of  UkUk  also  equals  k. 

The  SVD  solution  also  provides  a  formal  estimate  of  the  uncertainty  in  the 
solution  x  due  to  random  errors  in  the  data  (Wunsch,  1978).  The  relative  variance  of  the 
j'th  unknown  is  given  by 
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*    v2 
1-1   X2 


where  <Xj2>  is  the  solution  variance,  vjf  is  the  element  in  the  ift  column  and  y^  row  of  the 
V  matrix,  and  a2  is  the  a  priori  problem  variance.  Note  that  this  illustrates  the  trade-off 
between  resolution  and  variance,  in  that  as  k  is  increased,  smaller  and  smaller  values  of 
X,  are  kept  in  the  summation,  the  unknowns  become  more  well  resolved  but  the  solution 
variance  increases. 

4.4   Inverse  Results  for  Cruise  OC216 

The  inverse  model  for  the  December  cruise,  OC216  is  simpler  than  the  one  used 
for  the  April  cruise,  and  will  be  presented  first.  Two  inverse  calculations  were  made,  one 
to  balance  the  flow  between  legs  1  and  2  and  the  other  to  balance  the  flow  between  legs 
3  and  4.  In  chapter  3  we  noted  a  large  (7  Sv)  imbalance  between  legs  2  and  3  which  was 
deemed  to  be  due  to  edge  effects  caused  by  the  large  cold-core  ring  on  the  track  near 
Bermuda.  The  magnitude  of  this  change  in  transport  is  similar  to  that  reported  by  Halkin 
and  Rossby,  who  saw  shifts  in  Gulf  Stream  transport  of  order  10  Sv  in  a  7  day  period. 
Because  we  believe  this  change  in  transport  was  due  to  short-term  variability,  no  attempt 
was  made  to  balance  the  four  leg  transports  together. 

The  error  terms  used  in  the  row  weighting  for  cruise  OC205  were  estimated  to  be 
5  cm/s  for  the  individual  reference  velocities,  .2  Sv  for  the  overall  mass  transport,  1  Sv 
for  the  lower  4  layer  transports,  3  Sv  for  the  top  layer,  and  1.5  Sv  for  the  second  layer. 
The  shallower  layers  were  assigned  larger  errors  since  ageostrophic  velocity  components 
are  largest  near  the  surface. 
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Figure  4.1.  (a)  Plot  of  residual  norm  vs.  solution  norm  for  determining  solution  rank  for  legs  1  and  2  of  cruise 
OC216.   (b)  Plot  of  eigenvalues  vs.  rank  for  the  same  problem. 
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Cross-isopycnal  velocities  cruise  OC216  legs  1,2 
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Figure  4.2.   Cross-isopycnal  velocity  w\  vs.  rank  for  cniise  OC216  legs  1  and  2. 

In  order  to  obtain  the  "best"  solution  using  the  SVD,  the  rank  of  the  system  must 
be  estimated.  The  desired  rank  is  one  that  minimizes  the  residual  norm  without  making 
the  values  of  the  individual  velocity  corrections  too  large.  One  method  of  choosing  the 
optimum  rank  is  the  Levenberg-Marqardt  stabilization  method  (Lawson  and  Hanson, 
1974),  which  plots  the  residual  norm  vs.  the  solution  norm  for  each  rank.  The  SVD 
representation  of  the  Levenberg-Marquardt  diagram  for  legs  1  and  2  of  cruise  OC216  is 
shown  in  Figure  4.1  along  with  a  plot  of  the  eigenvalues  as  a  function  of  rank.  The 
optimum  rank  is  the  point  above  which  the  solution  norm  increases  sharply  without  a 
corresponding  drop  in  residuals.  This  should  correspond  to  a  drop  in  the  magnitude  of 
the  eigenvalues.  The  rank  of  this  problem  appears  to  be  approximately  97  or  98.  In 
order  to  further  evaluate  the  effect  of  rank  on  the  solution,  the  cross-isopycnal  velocities, 
w\  were  plotted  vs.  rank.  Figure  4.2    shows  that  the  magnitude  of  the  cross-isopycnal 
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Figure  4.3.     (a)  Plot  of  residual  norm  vs.  solution  norm    for  legs  3  and  4  of  cruise  OC216.     (b)  Plot  of 
eigenvalues  vs.  rank  for  the  same  problem. 
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velocities  increases  dramatically  above  rank  97.  Rank  97  was  chosen  as  the  optimum  rank 
of  the  system  based  on  keeping  the  cross-isopycnal  flow  small  while  balancing  the  layer 
and  total  transports  as  closely  as  possible.  A  similar  analysis  of  the  inverse  model  for 
legs  3  and  4  yields  an  optimum  rank  of  90   (Figure  4.3  and  Figure  4.4). 


Cross-lsopycnal  velocities  cruise  OC216  legs  3,4 
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Figure  4.4.   Cross-isopycnal  velocity  w\  vs.  rank  for  cruise  OC216  legs  3  and  4. 

The  velocity  corrections  for  the  rank  97  solution  for  legs  1  and  2  and  the  rank  90 
solution  for  legs  3  and  4  are  shown  in  Figure  4.5  with  error  bars.  As  discussed  by  Joyce 
et  al.  (1986),  the  error  bars  have  two  contributors:  (1)  that  due  to  noise  in  the  data  and 
(2)  that  due  to  failure  to  resolve  the  individual  value  of  q{  or  w*.  In  order  to  determine 
the  resolution  of  the  unknowns,  we  must  look  at  the  diagonal  elements  of  the  resolution 
matrix,  Vkv£.  These  values  are  plotted  vs.  unknown  in  Figure  4.6.  For  the  leg  1  and  2 
model,  the  first  100  unknowns  are  the  individual  reference  velocity  corrections,  qit  and 
the  last  six  unknowns  are  the  cross-isopycnal  velocities  w\.  We  can  see  that  most  of  the 
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Figure  4.5.  Reference  velocity  corrections  for  cruise  OC216.  Distances  are  measured  from  the  200  m  isobath. 
Some  representative  error  bars  are  provided.  The  error  bars  have  two  components,  data  noise  and  failure  to 
resolve.  Note  that  the  velocities  which  received  significant  corrections  are  well  resolved.  The  poorly  resolved 
velocities  (large  error  bars)  have    a  correction  of  zero. 


qK  are  well  resolved,  but  five  of  the  velocities  are  poorly  resolved.  Four  of  the  five  w\  are 
poorly  resolved.  The  velocity  corrections  required  to  balance  the  transports  are  much 
smaller  than  the  a  priori  estimate  of  ±5  cm/s  would  suggest.  For  each  leg,  significant 
adjustments  were  required  for  only  two  or  three  of  the  reference  velocities.  The  largest 
error  bars  in  Figure  4.5  correspond  to  the  velocities  which  have  not  been  resolved  by  the 
inverse  model.  The  rank  90  solution  for  legs  3  and  4  yields  similar  results.  Again,  most 
of  the  q{  are  well  resolved,  with  the  w\  in  general  being  less  well  resolved. 


Resolution  of  unknowns  OC216  legs  1  and  2 


Resolution  of  unknowns  OC216  legs  3  and  4 


Figure  4.6.  Diagonal  elements  of  VkV[  for  legs  1  and  2  rank  97  (top),  legs  3  and  4  rank  90  (bottom). 

Examining  the  diagonal  elements  of  the  UkUj  matrices,  we  find  that  for  the  leg 
1  and  2  model,  the  100  equations  for  the  reference  velocities  contribute  90.15  of  the  total 
rank  (97).  The  6  layer  transport  balance  equations  contribute  5.87,  and  the  total  transport 
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Cross-Isopycnal  velocities  between  OC216  legs  1  and  2 


Interface 


Figure  4.7.   Cross-isopycnal  velocities. 

equation  contributes  .98.  This  means  that  the  layer  and  total  mass  transport  balance 
equations  are  contributing  nearly  their  maximum  allowable  contributions  to  the  model, 
and  that  some  of  the  reference  velocity  equations  are  not  contributing.  The  results  of  the 
model  for  legs  3  and  4  are  similar,  with  the  92  equations  for  the  qi  contributing  83.1,  the 
6  layer  balance  equations  contributing  5.91,  and  the  total  mass  balance  equation 
contributing  .99  to  the  total  rank  of  90. 

The  values  of  w\  for  the  five  interface  areas  between  legs  1  and  2  and  between 
legs  3  and  4  are  plotted  in  Figure  4.7.  Due  to  the  small  values  and  the  lack  of 
resolution,  the  error  bars  are  large  enough  so  that  the  values  of  w\  are  not  statistically 
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distinguishable  from  zero.   The  transport  between  layers  due  to  these  velocities  is  very 
small  compared  to  the  horizontal  transports  across  the  sections. 

The  resulting  individual  layer  and  total  transports  are  given  in  Table  4.1.  The 
difference  between  these  values  and  the  values  given  in  chapter  3  are  small,  but  the 
formal  error  in  the  transport  calculations  has  been  reduced  from  -5  Sverdrups  to  less  than 
-1  Sverdrup  for  the  transports.  The  total  transports  for  each  leg  are  54.45  ±  .81  Sv  for 
leg  1,  54.67  ±  .32  Sv  for  leg  2,  63.07  ±  .40  Sv  for  leg  3,  and  62.39  ±  .62  Sv  for  leg  4. 

Plots  of  accumulated  transport  as  a  function  of  distance  from  the  200  meter 
isobath  are  given  in  Figure  4.8,  comparing  the  results  of  the  inverse  model  with  the 
transports  computed  from  the  direct  addition  of  reference  velocity  to  the  geostrophic 
profiles  in  chapter  3.  Due  to  the  small  changes  in  the  reference  velocities,  the  changes 
in  accumulated  transport  across  the  sections  are  quite  small,  but  were  enough  to  balance 
the  overall  transport  between  section  1  and  2  and  between  sections  3  and  4  to  within  1 
Sverdrup. 
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Table  4.1.   Transports  for  cruise  OC216  from  inverse  model. 


OC216  Layer  and  Total  Section  Transports 


Layer 

Leg  1 

Leg  2 

Leg  3 

Leg  4 

1 

22.449 

22.420 

22.132 

26.531 

2 

9.554 

10.869 

14.355 

12.864 

3 

8.572 

9.027 

11.495 

9.096 

4 

7.033 

6.816 

8.880 

7.717 

5 

5.397 

5.063 

5.376 

4.798 

6 

5.063 

0.476 

0.833 

1.388 

Total 

54.451 

54.670 

63.071 

62.393 

Rank  97  inverse  solution  for  legs  1,2. 
Rank  90  inverse  solution  for  legs  3,4. 


Cross-isopycnal  Transport 

Interface 

w\  12 

Tran  12 

w]  34 

Tran  34 

1 

0.198 

0.06 

0.882 

0.28 

2 

0.032 

0.01 

0.343 

0.11 

3 

-0.005 

-0.002 

-0.251 

-0.08 

4 

0.056 

0.005 

0.059 

0.006 

5 

-0.412 

-0.03 

-0.757 

-0.06 

w\  in  m/s  x  10"5 

Transport  in  Sverdrups  (kg/s  x  109) 

Cross-isopycnal  transports  computed  by  multiplying  the  cross-isopycnal  velocity  w\ 

by  the  area  of  the  interface. 
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Accumulated  transport  OC216  leg  1  (top),  leg  2  (bottom) 
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Figure  4.8.  Transport  as  a  function  of  distance  for  the  four  sections  of  cruise  OC216.  The  solid  lines  are 
transports  from  direct  combination  of  ADGP  and  geostrophic  data.  The  dotted  lines  are  the  result  of  the  inverse 
models. 
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4.5   Inverse  Results  for  Cruise  OC205 

The  geometry  of  the  April  cruise  allows  us  to  balance  the  mass  transport  between 
three  legs.  For  the  center  track  between  Bermuda  and  Cape  Cod,  the  second  transit  (leg 
4)  was  chosen  to  balance  with  legs  3  and  6,  since  it  is  the  leg  closest  in  time  to  legs  3 
and  6.  In  this  model  there  is  a  set  of  equations  balancing  the  flow  in  leg  3  with  the  flow 
in  leg  4,  and  another  set  of  equations  relating  the  flow  in  leg  4  to  the  flow  in  leg  6.  The 
unknowns  in  this  case  are  the  reference  velocities  q,  for  the  three  legs,  and  the  cross- 
isopycnal  velocities  w\  for  the  areas  between  legs  3  and  4  and  between  legs  4  and  6.  The 
number  of  equations  for  this  model  is  166,  with  152  equations  for  the  reference  velocities, 
12  individual  layer  balance  equations  (6  for  the  leg  3-6  balance,  and  6  for  the  leg  4-6 
balance),  and  2  total  section  mass  balance  equations.  The  number  of  unknowns  is  162 
(152  reference  velocities  +  10  cross-isopycnal  velocities).  The  initial  mass  transports  for 
this  model  are  significantly  more  out  of  balance  than  the  ones  for  the  OC216  cruise.  The 
row  weighting  was  changed  for  this  model,  allowing  an  error  in  total  transport  of  ±  7  Sv, 
while  holding  the  individual  transports  to  ±3  Sv  for  the  top  layer,  ±1.5  Sv  for  the  second 
layer,  and  ±1  Sv  for  the  bottom  4  layers.  As  we  will  see  when  examining  the  UkUj 
elements,  this  reduces  the  weight  of  the  total  transport  equations.  By  evaluating  the 
Levenberg-Marquardt  diagram,  the  plot  of  eigenvalues  vs.  rank,  and  the  plot  of  cross- 
isopycnal  velocity  vs.  rank,  (Figure  4.9  and  Figure  4.10),  the  optimum  rank  for  the  SVD 
of  this  model  was  chosen  to  be  152. 

Examining  the  resolution  matrix  (Figure  4.11),  we  find  that  the  reference  level 
velocities  qt  are  nearly  fully  resolved,  while  the  w\  are  essentially  unresolved.  The 
velocity  corrections  (Figure  4.12)  are  much  larger  than  those  seen  in  the  OC216  models. 
In  order  to  balance  the  transports,  the  model  has  added  additional  bias  to  each  leg,  as  well 
as  changing  individual  velocities.  The  the  mean  values  of  the  velocity  corrections  in 
Figure  4.12  are  -1.23  cm/s  for  leg  3,  -.32  cm/s  for  leg  4,  and  +.85  cm/s  for  leg  6. 
Examining  the  values  of  w*  (Figure  4.13),  we  find  that  since  the  cross-isopycnal  velocities 
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Norms  for  SVD  of  OC205  legs  3,  4,  and  6 
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Figure  4.9.     (a)  Levenberg-Marquardt  diagram  for  3  leg  inversion  model  for  cruise  OC205.     (b)  Plot  of 
eigenvalue  magnitude  vs.  rank  for  the  same  model. 
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Cross-lsopycnal  velocities  cruise  OC205  legs  3,4 
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Figure  4.10.   Cross-isopycnal  velocities  vs.  rank. 
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Figure  4.11.  Diagonal  elements  of  resolution  matrix  VkVj  for  Rank  152. 

have  not  been  well  resolved,  the  error  bars  are  much  larger  than  the  velocities  themselves. 
The  w\  are  therefore  indistinguishable  from  zero. 

The  diagonal  elements  of  UkUj  show  that  the  contribution  of  the  equations  to  the 
solution  is  different  from  the  OC216  model.  The  152  reference  velocity  equations 
contribute  147.94,  the  12  equations  for  the  layer  transports  (six  for  legs  3-4,  six  for  legs 
4-6)  contribute  3.60,  and  the  two  total  mass  balance  equations  contribute  .46  to  the  rank 
of  152.  In  this  case,  little  weight  was  put  on  the  total  mass  balance  equations  due  to  the 
large  errors  allowed  for  these  equations.  The  individual  layer  equations  also  contribute 
only  fractionally  to  the  solution,  with  most  of  the  rank  coming  from  the  reference  velocity 
equations. 

The  individual  layer  and  total  transports  for  the  inverse  rank  152  solution  are 
summarized  in  Table  4.2.  Despite  the  low  weight  on  the  total  transport  constraint,  the 
model  has  lowered  the  total  transport  differences  to  4  Sv.  between  legs  3  and  4  and  2  Sv. 
between  legs  4  and  6.  The  total  transports  are  now  54.56  ±1.3  Sv.  for  leg  3,  50.7  ±.87 
Sv.  for  leg  4,  and  49.69  ±  1 .29  Sv.  for  leg  6.  The  cross-isopycnal  velocities  are  of  order 
10"5  m/s,  which  is  an  order  of  magnitude  higher  than  seen  in  the  OC216  solutions,  but 
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Figure  4.12.  Velocity  corrections  for  cruise  OC205  Rank  152.  Error  bars  include  error  due  to  noise,  and  error 
due  to  failure  to  fully  resolve. 

still  small  enough  that  cross-isopycnal  velocity  causes  only  small  changes  of  order  .5  Sv 
in  layer  transport  between  legs. 
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Figure  4.13.   Cross- isopycnal  velocities  cruise  OC205  rank  152. 

The  transports  calculated  from  the  inverse  model  are  shown  compared  with  the 
transports  calculated  in  chapter  3  in  Figure  4.14.  The  changes  in  transport  introduced  by 
the  inverse  solution  are  significantly  larger  than  those  produced  for  the  OC216  cruise, 
since  the  initial  transports  had  much  larger  differences  than  the  inputs  to  the  OC216 
model.  Leg  4  required  the  smallest  amount  of  change,  which  is  not  surprising,  since  its 
transport  was  initially  between  the  other  two  legs.  The  transport  for  leg  3  has  been 
reduced,  and  the  transport  for  leg  6  increased  to  more  nearly  match  leg  4.  The  variability 
of  the  transport  near  Bermuda  can  be  seen  by  plotting  the  transports  of  the  three  legs  vs. 
distance  from  the  Bermuda  vertex  (Figure  4.15).  The  tracks  diverge  with  distance  from 
Bermuda,  but  might  all  be  expected  to  have  the  same  type  of  transport  change  with 
distance  near  the  point  where  the  tracks  meet.  We  can  see  that  this  is  not  the  case, 
indicating  the  currents  near  Bermuda  have  changed  during  the  time  between  legs. 
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Table  4.2.    Transports  for  cruise  OC205  legs  3,  4,  and  6  from  inverse  model. 


OC205  Layer  and  Total  Section  Transports 


Layer 

Leg  3 

Leg  4 

Leg  6 

1 

22.36 

17.48 

16.53 

2 

12.70 

10.02 

7.93 

3 

8.62 

10.40 

10.04 

4 

8.08 

8.19 

10.34 

5 

2.49 

5.35 

6.19 

6 

0.31 

-0.74 

-1.33 

Total 

54.56 

50.70 

49.69 

Inverse  model  rank  152 

Transports  in  Sverdrups  (kg/s  x  10"9) 


Cross-Isopycnal  Transport 


Interface 

w*  34 

Tran  34 

w*   46 

Tran  46 

1 

1.62 

0.59 

-1.10 

-0.21 

2 

-1.51 

-0.55 

0.07 

0.01 

3 

-0.42 

-0.15 

2.80 

0.54 

4 

-2.54 

-0.50 

-0.27 

-0.04 

5 

1.69 

0.20 

-2.02 

-0.21 

w*  in  m/s  x  10*3 

Transports  in  Sverdrups  (kg/s  x  10"9) 

Cross-isopycnal  transports  computed  by  multiplying  w\  by  the  interface  area. 
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Figure  4.14.  Transport  as  a  function  of  distance  for  OC205  legs  3,  4,  and  6.  The  solid  lines  are  transport  from 
direct  combination  of  ADCP  and  geostrophic  data  from  chapter  3.  The  dotted  lines  are  the  result  of  the  inverse 
model. 
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Accumulated  Transports  OC205  legs  3,  4,  and  6 

Inverse  solution  rank  152 
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Figure  4.15.  Transports  OC205  legs  3,  4,  and  6  vs.  distance  from  Bermuda. 

This  graph  also  illustrates  the  change  in  the  Gulf  Stream  transport  as  it  moves 
downstream.  The  sharp  rise  in  transport  from  the  Gulf  Stream  provides  approximately 
55  Sv.  in  the  upper  765  meters  for  both  legs  3  and  4,  but  has  increased  to  65  Sv.  by  leg 
6.  These  estimates  are  very  sensitive  to  the  definition  of  what  portion  of  the  sections 
constitutes  the  "Gulf  Stream". 
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Chapter  5 


Temporal  Variability  between  Bermuda  and  Cape  Cod 


5.1    Gulf  Stream  Variability 


The  Acoustic  Doppler  Current  Profiler  (ADCP)  allows  us  to  compare  velocity 
profiles  computed  using  the  dynamic  method  to  directly  measured  velocity  profiles.  Since 
it  measures  both  horizontal  components  of  velocity,  we  can  also  directly  calculate  several 
of  the  terms  in  the  horizontal  momentum  equations  which  are  neglected  when  using  the 
geostrophic  approximation. 

The  Altimeter  track  between  Bermuda  and  Cape  Cod  was  traversed  three  times 
during  the  OC205  cruise,  and  twice  during  the  OC216  cruise.  By  examining  the  ADCP 
velocities  for  the  different  legs  along  with  AVHRR  sea  surface  temperature  images,  we 
can  see  how  the  Gulf  Stream  is  changing  during  the  cruises.  During  the  April  cruise, 
OC205,  the  Gulf  Stream  moved  to  the  south  and  underwent  a  rotation  in  the  clockwise 
direction  during  the  cruise,  due  to  the  trough  of  a  meander  which  crossed  the  track  during 
the  cruise.  This  rotation  is  seen  in  the  100  meter  doppler  vectors  for  the  three  tracks 
(Figure  5.1).  The  position  of  the  meander  relative  to  the  track  on  various  days  can  be 
seen  in  the  AVHRR  images  (Figure  5.2).  During  the  December  cruise  (OC216),  a  small 
southward  shift  in  the  position  of  the  Gulf  Stream  was  observed  between  legs  1  and  4, 
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Figure  5.1.  ADCP  velocities  at  100  meters  depth,  cruise  OC205  legs  1,  4,  and  7.  Leg  1  has  been  moved  2° 
west,  leg  7  has  been  moved  2°  east 

and  additional  variability  was  observed  from  a  large  cold-core  ring  which  moved  across 
the  track  (Figure  5.3).  Unfortunately,  due  to  heavy  cloud  cover,  no  AVHRR  images  of 
the  sea  surface  are  available  for  the  December  cruise. 
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Figure  5.2.    Infrared  AVHRR  images  for  days  94  (upper  left),  102  (upper  right),  108  (lower  left),  and  111 
(lower  right),  with  GEOSAT  track. 
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Figure  5.3.    ADCP  velocities  at  1(X)  meters  depth  for  cruise  OC216  legs  1  and  4.   Leg  4  has  been  shifted  2° 
east. 

5.2    Ageostrophic  effects 

The  multiple  crossings  allow  us  to  evaluate  time  as  well  as  space  derivatives  of 
velocity  along  the  Cape  Cod-Bermuda  track.  Velocities  and  their  derivatives  were 
measured  in  a  rotated  coordinate  system,  where  x  is  the  across-track  direction,  and  y  is 
the  along  track  direction.   Velocity  in  the  x  direction  is  u  (positive  to  the  Northeast)  and 
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velocity  in  the  y  direction  is  v  (positive  to  the  Northwest).    The  y  direction  (ship's 
heading)  is  334°. 

The  horizontal  momentum  equations  in  the  rotated  coordinate  system  are: 

ut  +  uux  +  vu    +  wuz  -  fv    =    -Px  (5.1) 

v    +  uv    +  vv    +  wv    +  fu     =     -P  (5.2) 

x  x  y  z  y 

In  order  to  examine  the  terms  in  the  equations  at  the  same  points  in  space,  the 
ADCP  vectors  were  filtered  onto  a  straight  line  track  using  a  Gaussian  filter  with  a  filter 
half-width  of  10  km.  A  filtered  velocity  vector  was  generated  every  10.2  km  between 
32°-30'  N  and  40°  N.  Local  time  derivatives  were  generated  for  each  point  by  taking  the 
difference  in  u  and  v  between  legs  and  dividing  by  the  time  between  measurements. 
Spatial  derivatives  were  obtained  at  the  same  points  as  the  time  derivatives  using  a 
centered  difference  scheme. 

Spatial  derivatives  are  only  directly  measurable  in  the  along  track  direction.  A 
quasi-geostrophic  approximation  (c.f.  Pedlosky,  1979)  can  be  used  to  justify  neglect  of 
the  vertical  velocity  in  the  advection  of  horizontal  momentum  and  in  the  divergence 
equation.  However,  in  the  Gulf  Stream,  vertical  velocities  can  be  quite  large.  Bower 
(1989)  has  shown  that  vertical  velocity  can  be  related  to  vorticity  changes  following  a 
"Lagrangian"  RAFOS  float.   Using  the  continuity  equation,  it  is  possible  to  estimate  ux. 

This  approximation  is  justified  (Joyce,  Kelly,  Schubert,  and  Caruso,  1990)  using  values 
for  wz  in  the  Gulf  Stream  of  10"*  s"1  from  Bower  and  using  measured  values  for  v  . 

[wj/[vy]  =  lO^/lO5  =  0.1 
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The  vertical  advection  term  can  also  be  shown  to  be  small,  using  values  for  w  of  5  x  10"4 
m/s  from  Bower  and  measured  values  for  the  remaining  terms. 


[wuj/[vuy]  =  2  x  10^/10°  =  0.2 

Using  the  above  scaling,  equation  (1)  can  be  rewritten  as 

u.  +  vu    -  uv    -  fv     =     -P,  -  w«,  +  uw,  (5.4) 

t  y  y       J  x  z  z 

« 

where  all  the  terms  on  the  left-hand  side  can  be  directly  measured,  and  the  terms 
involving  w  on  the  right-hand  side  are  small. 

During  the  last  leg  of  cruise  OC216,  a  short  section  was  made  perpendicular  to 
the  ship  track.  This  allowed  measurement  of  derivatives  in  the  x  direction,  providing  an 
opportunity  to  check  the  above  scaling  arguments  at  one  point  in  the  Gulf  Stream.  The 
cross-track  was  approximately  90  km  long,  and  crossed  the  satellite  track  about  five  km 
south  of  the  peak  velocity  point  of  the  satellite  track.  At  the  point  where  the  two  tracks 
cross,  a  comparison  can  be  made  of  ux  and  vy.  Using  the  same  Gaussian  filter  described 
previously,  the  three  values  of  ux  and  vy  closest  to  the  crossing  point  have  been  averaged 
at  50,  100,  150,  and  200  meters.  The  value  of  vvz  at  each  of  these  depths  has  been 
computed  using  the  continuity  equation.  These  values  with  standard  deviations  are 
tabulated  in  Table  5.1. 
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Table  5.1.   Velocity  derivatives  in  the  continuity  equation. 


Depth 

u, 

a* 

vy 

°y 

w* 

Vz 

50 

-1.239 

2.442 

-5.015 

1.568 

6.254 

2.902 

100 

0.193 

1.750 

-1.031 

1.790 

0.838 

2.503 

150 

-1.734 

2.819 

0.005 

1.642 

1.729 

3.262 

200 

0.021 

3.629 

-1.239 

2.918 

1.260 

4.657 

Depths  are  in  meters,  derivatives  in  s"1  x  10"6. 

Standard  deviations  in  x  and  y  directions  are  of  the  three  data  points  closest  to  the  point  where  the  tracks 
cross. 

a,-(o,J  +  0yT 


100  150 

Depth    (meters) 

Figure  5.4.   Velocity  derivatives  du/dx  (*)  and  -dv/dy  (o)  at  various  depths,   (sec"1  x  10"3) 
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Values  of  u,  and  -vy  are  plotted  with  error  bars  in  Figure  5.4.  For  the  100  to  200 
meter  values,  the  difference  between  the  values  of  ux  and  -vy  falls  within  the  error  bars. 
There  is  a  large  difference  between  the  terms  at  50  meters,  indicating  that  wz  is 
measurably  different  from  zero  at  shallow  depths.  This  agrees  with  Bower's  finding  that 
wz  tends  to  have  larger  values  at  shallower  density  surfaces. 
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Figure  5.5.   NOAA  Oceanographic  Analysis  Chart  dated  11  Dec  1989.  The  ship's  tracks  have  been  added  as 
heavy  dashed  lines.   The  Gulf  Stream  is  indicated  by  "GS\ 

The  positive  value  of  wz  at  50  meters  indicates  a  horizontal  divergence,  which 
Bower  found  to  occur  in  the  right  side  of  the  stream  when  downstream  of  a  trough,  and 
on  the  left  side  of  the  stream  when  downstream  of  a  meander  crest.  Our  perpendicular 
section  crossed  the  track  to  the  south  of  the  peak  in  Gulf  Stream  velocity,  indicating  that 
our  data  point  for  wz  is  in  the  right  side  of  the  stream,  and  therefore  are  downstream  of 
a  trough.   There  are  no  AVHRR  images  available  to  confirm  the  orientation  of  the  Gulf 
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Stream,  however  the  weekly  NOAA/National  Ocean  Service  weekly  Oceanographic 
analysis  chart  dated  1 1  Dec  1989  (Figure  5.5)  indicates  that  the  Gulf  Stream  upstream  of 
the  ship's  track  is  relatively  straight,  unlike  the  situation  in  the  April  cruise. 

Using  the  value  of  wz  at  50  meters  and  assuming  that  w=0  at  the  surface,  we  can 
estimate  the  value  of  w  at  50  meters  in  the  Gulf  Stream. 

w50    =    -^oz    =    (6.254xl0"6)(-50)    =     -3.1X10"4  mfs  (5.5) 

etc 


Actual  measured  values  of  w  in  the  Gulf  Stream  are  of  order  8  x  10"*  m/s  rms 
(Bower  and  Rossby,  1989)  at  depths  between  300  and  700  meters.  Our  computed 
velocity  is  smaller  than  these  values,  but  is  measured  at  a  much  shallower  depth. 
If  the  sign  of  wz  remains  positive  (although  decreasing)  with  depth,  the  magnitude  of  w 
would  increase  with  depth.  Hall  (1986)  estimates  this  depth  of  maximum  w  to  be  875 
meters,  based  on  estimates  of  vertical  velocity  at  a  moored  current  meter  site.  The  above 
data  do  not  support  approximating  ux  by  -vy,  since  vvz  is  of  the  same  order  of  magnitude 
as  the  other  terms.  It  is  accurate  to  say  that  ux  will  be  of  the  same  order  of  magnitude 
as  vy  when  deciding  which  terms  to  neglect  in  the  momentum  equations. 

A  clear  picture  of  the  relative  magnitudes  of  the  terms  of  the  momentum  equations 
can  be  made  by  simply  graphing  all  of  the  measurable  terms.  Figure  5.6  contains  graphs 
of  the  momentum  equation  terms  at  100  and  200  meters  for  cruise  OC205.  The 
momentum  equation  terms  for  cruise  OC216  at  50,  100,  150,  and  200  meters  are  shown 
in  Figure  5.7.  In  each  plot,  the  terms  are  average  values  for  the  crossings  conducted  on 
that  cruise;  three  crossings  for  the  first  cruise,  and  two  for  the  second  cruise.  The  curves 
change  little  with  depth,  indicating  that  uz  and  vr  are  small. 
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X  momentum  equation  terms  at  100  meters 
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Figure  5.6.   Momentum  equation  terms  for  cruise  OC205. 
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X  momentum  equation  terms  at  200  meters 
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Y  momentum  equation  terms  at  200  meters 
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X  momentum  equation  terms  at  50  meters 
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Figure  5.7.   Momentum  equation  terms  for  cruise  OC216. 
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X  momentum  equation  terms  at  100  meters 
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X  momentum  equation  terms  at  150  meters 
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Y  momentum  equation  terms  at  150  meters 
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X  momentum  equation  terms  at  200  meters 


Y  momentum  equation  terms  at  200  meters 
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Looking  at  the  relative  magnitudes  of  the  terms  in  the  x  equations,  we  see  that 
the  local  acceleration  term  wt  is  the  smallest  with  magnitude  =10"*  m/s2.  The  convective 
acceleration  terms  uvy  and  vuy  have  maximum  values  of  1-2  x  10"5  m/s2,  and  the  Coriolis 
term  -fv  has  a  maximum  value  of  5  x  10"5  m/s2.  Using  the  geostrophic  approximation  in 
the  cross-track  direction  would  give  velocity  errors  of  20  to  35  percent. 

Adding  all  the  measured  ageostrophic  terms  together,  we  can  get  a  feel  for  what 
the  pressure  gradient  must  be.  The  sums  of  x  momentum  equation  terms  ut  -  uvy  +  vuy 
and  y  momentum  equation  terms  vt  +  vvy  with  error  bars  are  shown  in  Figure  5.8  for 
cruise  OC205  and  in  Figure  5.9  for  cruise  OC216.  For  cruise  OC205,  the  sum  of  «,  -  uvy 
+  vuy  in  the  Gulf  Stream  (latitude  37°  -  38°  N  in  Figure  5.8)  is  positive,  indicating  that 
the  pressure  gradient  in  the  x  direction  would  have  to  be  be  larger  than  required  to 
balance  the  Coriolis  term.  This  additional  pressure  gradient,  if  not  attenuated  with  depth 
would  tend  to  push  water  to  the  south  at  speeds  of  approximately  10  cm/s  at  depths  where 
the  non-linear  terms  decay,  because  the  Coriolis  term  is  not  large  enough  in  itself  to 
balance  the  pressure  gradient.  This  agrees  with  the  findings  of  Bower  and  Rossby  (1989) 
that  cross-frontal  velocity  would  be  to  the  south  downstream  of  the  crest  of  a  Gulf  Stream 
meander. 

Looking  at  the  sum  of  the  x  momentum  equation  terms  for  cruise  OC216,  the 
results  are  less  clear  in  the  Gulf  Stream.  There  is  a  large  negative  spike  at  38.5°  N,  right 
in  the  center  of  the  Stream.  The  overall  sense  of  the  ageostrophic  terms  is  small,  which 
means  the  overall  pressure  gradient  across  the  stream  should  be  about  the  same  as  that 
required  to  balance  the  Coriolis  term.  Since  it  appears  that  on  this  cruise  we  crossed  a 
fairly  flat  section  of  the  stream,  with  no  meander,  this  is  consistent  with  Bower  and 
Rossby. 

The  graphs  of  the  y  momentum  equation  terms  show  only  three  terms,  v„  vvy,  and 
fit.   The  term  which  we  are  unable  to  measure,  mvx  is  known  as  the  cyclostrophic  term, 
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Sum  of  x  equation  terms  at  100  meters  Cruise  OC205 
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Sum  of  y  equation  terms  at  100  meters  Cruise  OC205 


Figure  5.8.  Sums  of  momentum  equation  terms  at  100  meters,  cruise  OC205.  Error  bars  are  standard  deviation 
of  u,  -  uv  +  vu^  for  the  x  equation,  and  standard  deviation  of  v,  +  wy  for  the  y  equation. 
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X  momentum  equation  terms  at  100  meters  Cruise  OC216 


6 

X 


4h 

3 

2 

1 

0 
-1 
-2 
-3 
-4 


-5 


40 


39 


38 


34 


37  36  35 

Latitude 
Y  momentum  equation  terms  at  100  meters  Cruise  OC216 


33 


32 


Figure  5.9.  Sums  of  momentum  equation  terms  at  100  m  depth,  cruise  OC216.  Error  bars  are  standard  deviation 
of  u,  -  wv  +  viiy  for  the  x  equation,  and  standard  deviation  of  vt  +  wy  for  the  y  equation. 
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which  can  be  large  if  there  is  significant  curvature  in  the  Gulf  Stream,  which  would  cause 
vx  to  be  significant.  From  the  AVHRR  pictures  during  cruise  OC205,  we  know  that  we 
are  in  a  Gulf  Stream  meander,  and  can  estimate  uvx  using  the  radius  of  curvature.  For 
the  curvature  present  during  cruise  OC205,  the  assumption  of  pure  geostrophy  causes  u 
to  be  overestimated  by  approximately  15%  in  the  center  of  the  Gulf  Stream.  A  more 
detailed  analysis  of  this  effect  for  cruise  OC205  is  conducted  in  the  next  chapter. 

For  the  OC216  cruise,  a  value  of  uvx  can  be  determined  for  one  point  in  the  Gulf 
Stream  using  the  data  from  the  cross-track  conducted  during  leg  4.  The  value  of  uvx  at 
this  point  is  -2  x  10"5  m/s2  at  a  depth  of  100  meters.  Adding  this  term  to  the  sum  of  v, 
+  vvy  shown  in  Figure  5.9,  the  sum  of  the  ageostrophic  terms  in  the  y  momentum 
equation  would  be  of  the  order  1-2  x  10"5  m/s2  near  the  center  of  the  Gulf  Stream  during 
cruise  OC216.  With  values  of  fu  in  the  Gulf  Stream  of  order  1.3  x  10"4  m/s2,  using 
geostrophy  to  compute  cross-track  velocity  u  would  result  in  errors  of  order  10%. 
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Chapter  6 

Comparison  of  Sea  Surface  Height  and  Velocity  to 

Satellite  Data 


The  GEOSAT  altimeter  measures  collinear  profiles  of  sea  surface  height  every 
seventeen  days.  By  calculating  dynamic  height  and  sea  surface  velocity  from  the 
shipboard  transits  along  one  of  the  GEOSAT  tracks,  we  are  able  to  get  an  independent 
estimate  of  sea  surface  height  to  compare  with  the  altimeter  data.  In  order  to  compare 
the  two  data  sets,  the  sea  surface  height  variations  due  to  the  geoid  must  first  be  removed 
from  the  altimeter  data.  A  mean  sea  surface  topography  is  achieved  by  averaging  many 
GEOSAT  passes  together.  Unfortunately,  subtracting  out  this  mean  signal  from  the  height 
data  also  subtracts  out  the  mean  dynamic  height.  Kelly  and  Gille  (1990)  developed  a 
method  to  add  a  synthetic  mean  Gulf  Stream  height  profile  back  into  the  altimeter 
residual  height  data  to  produce  a  total  dynamic  height  profile.  The  mean  Gulf  Stream  is 
modeled  with  a  Gaussian  profile  of  the  form 

(y-b)2 
u(y)    =    aexp\--x—f-    ,  (6.1) 

2c2 
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where  u  is  the  downstream  velocity  as  a  function  of  cross-stream  location  (y),  a  is  the 
peak  velocity,  b  is  the  location  of  the  center  of  the  stream,  and  c  is  a  width  parameter. 
The  dynamic  height  is  then  given  by  the  integral 

h(y)    =    -^fu(y)dy/.  (6.2) 

8  J 


The  mean  value  of  h(y)  when  added  to  the  residual  height  data  from  a  particular 
satellite  pass  will  give  an  estimate  of  the  sea  surface  dynamic  height  for  that  pass.  This 
height  can  be  differentiated  using  the  geostrophic  relation  to  get  a  sea  surface  velocity 
profile.  There  are  several  errors  introduced  by  using  this  method.  First,  the  GEOSAT 
tracks  are  not  exactly  collinear.  The  satellite  can  deviate  from  its  nominal  track  by  one 
kilometer  or  more  in  the  cross-stream  direction.  This  means  that  when  the  mean  sea 
surface  topography  is  calculated,  cross-stream  spatial  variations  in  sea  surface  height  will 
show  up  as  temporal  variations.  This  effect  is  most  pronounced  near  Bermuda,  where 
there  are  strong  cross-track  gradients  of  the  geoid. 

Another  restriction  of  the  model  is  that  mean  values  for  flows  outside  the  Gulf 
Stream  are  not  modeled.  The  mean  recirculation  north  and  south  of  the  Gulf  Stream  and 
any  mean  currents  associated  with  rings  are  not  taken  into  account. 

The  third  source  of  error  which  affects  the  computed  surface  velocity  comes  from 
assuming  geostrophy  at  the  sea  surface.  As  shown  in  the  previous  section,  there  are 
ageostrophic  components  in  the  momentum  equations  which  are  significant  in  the  surface 
layer.  In  regions  of  significant  curvature  in  the  Gulf  Stream,  the  "cyclostrophic"  term 
should  be  included  in  the  momentum  balance. 

Sea  surface  velocity  profiles  were  generated  from  each  of  the  three  shipboard 
tracks  from  cruise  OC205  and  for  the  closest  two  GEOSAT  passes.  The  ship's  track 
surface  velocities  are  from  the  combined  doppler  and  geostrophic  profiles.  Unfortunately, 
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Figure  6.1.    Sea  surface  velocities  and  sea  surface  heights  for  ship  and  satellite  tracks  between  Bermuda  and 
Cape  Cod. 


94 


none  of  the  shipboard  crossings  of  the  Gulf  Stream  coincided  in  time  with  a  satellite  pass. 
The  GEOSAT  pass  occurring  during  the  cruise  was  on  April  17,  between  legs  4  and  7  of 
the  ship's  track.  These  cross-track  sea  surface  velocity  profiles  are  presented  as  a  time 
series  in  Figure  6. 1 .  The  satellite  profiles  and  the  shipboard  profiles  show  a  consistent 
picture  of  the  Gulf  Stream  moving  southward  during  the  period  of  the  cruise,  with  the 
peak  cross-track  velocity  in  the  Gulf  Stream  decreasing  towards  the  end  of  the  cruise  due 
to  rotation  of  the  Gulf  Stream  axis.  The  values  of  peak  velocity  (about  2  m/s)  and  width 
of  the  Gulf  Stream  are  similar  in  both  satellite  and  shipboard  profiles,  although  the  peak 
Gulf  Stream  velocities  are  somewhat  higher  in  the  satellite  profiles.  The  Gulf  Stream 
appears  to  have  broadened  and  slowed  down  by  leg  7  because  the  meander  which  was 
crossed  on  this  leg  has  rotated  the  Gulf  Stream  axis  to  the  south.  The  GEOSAT  velocity 
profile  for  day  107  shows  the  broadening  associated  with  this  rotation,  but  not  the  dip  in 
peak  velocity.  Part  of  this  inconsistency  may  be  due  to  omission  of  the  cyclostrophic 
term  in  the  momentum  balance. 

The  AVHRR  infrared  satellite  images  presented  in  chapter  5  show  the  evolution 
of  the  Gulf  Stream  meander  which  was  located  on  the  satellite  track  at  the  time  of  the 
cruise.  The  meander  moved  from  east  to  west  and  the  trough  of  the  meander  deepened 
during  the  cruise.  Due  to  the  significant  curvature  in  the  Gulf  Stream  induced  by  this 
meander,  the  cyclostrophic  term  is  significant  in  the  momentum  balance.  The  momentum 
balance,  including  the  cyclostrophic  term  is 

t  +  *  -  -*£  •  (<UI 

R  *» 


where  u  is  the  velocity  in  the  direction  tangent  to  the  radius  of  curvature,  and  n  is 
distance  normal  to  the  radius  of  curvature.  In  the  track  oriented  coordinate  system  the 
equation  becomes: 
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RcosQ 


fu    =     -g 


dy 


■     fU8> 


(6.4) 


where  ur  is  the  cross-track  velocity  in  the  rotated  system,  0  is  the  angle  between  the 
downstream  direction  and  the  cross-track  direction,  and  wg  is  the  velocity  which  would 
exist  if  the  height  gradient  were  due  entirely  to  geostrophy. 
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Figure  6.2.  Gulf  Stream  curvature  at  times  of  day  90  and  day  107  GEOSAT  passes. 

At  the  time  of  the  day  90  satellite  pass,  the  radius  of  curvature  was  =170  km  and 
the  angle  between  the  Gulf  Stream  direction  and  the  normal  to  the  satellite  track  was  3°. 
Using  equation  (6.4),  the  geostrophic  peak  velocity  in  the  Gulf  Stream  of  2.2  m/s 
corresponds  to  a  peak  velocity  of  1.95  m/s  when  the  curvature  is  taken  into  effect.  For 
the  second  satellite  pass,  using  R  =  130  km,  6=  33°,  and  v^  =  2.05  m/s,  the 
corresponding  peak  velocity  including  the  cyclostrophic  term  is  1.74  m/s.    (Figure  6.2) 


96 


The  surface  velocity  computed  by  differentiating  the  sea  surface  slope  using 
geostrophy  is  about  .25  to  .3  m/s  too  high  at  the  peak  value  in  the  Gulf  Stream  given  the 
curvature  seen  at  the  time  of  the  satellite  passes.  Since  the  shipboard  velocities  are 
referenced  to  directly  measured  ADCP  velocities,  the  surface  velocities  computed  from 
shipboard  data  should  be  close  to  the  actual  velocity,  even  though  geostrophy  was  used 
to  compute  the  velocity  change  from  100  meters  to  the  surface. 

Smaller  scale  variations  in  velocity  also  appear  consistently  between  the  altimeter 
and  the  shipboard  profiles.  A  cold  core  ring  centered  at  35.6°  N  can  be  seen  in  both  the 
day  90  GEOSAT  and  the  leg  1  shipboard  profile.  The  large  fluctuations  in  the  altimeter 
velocities  south  of  33°  N  are  noise  due  to  geoid  spatial  variability  near  Bermuda  as 
discussed  earlier. 

More  insight  can  be  gained  by  examining  a  time  series  of  sea  surface  heights.  A 
time  series  of  height  profiles  from  the  three  ship  crossings  and  two  altimeter  crossings  are 
shown  in  the  lower  half  of  Figure  6.1  (page  99).  In  this  case,  the  heights  measured  by 
the  altimeter  are  the  direct  measurement,  and  the  velocities  measured  from  the  ship  are 
integrated  to  give  a  geostrophic  dynamic  height  profile.  The  height  differences  across  the 
Gulf  Stream  determined  geostrophic  ally  from  the  shipboard  data  are  low  due  to 
neglecting  the  cyclostrophic  term.  In  order  to  see  how  close  the  best  shipboard  estimate 
of  dynamic  height  would  come  to  the  altimeter  height,  the  velocity  data  from  legs  4  and 
7  were  interpolated  to  the  time  of  the  GEOSAT  pass.  A  cyclostrophic  correction  was 
determined  for  each  interpolated  doppler  vector  using  the  component  of  velocity  in  the 
radial  direction  and  the  distance  from  each  vector  position  to  the  center  of  the  radius  of 
curvature.  The  100  meter  doppler  vectors  were  used,  assuming  that  the  cyclostrophic 
term  would  not  be  significantly  different  at  100  meters  than  at  the  surface.  The 
cyclostrophic  correction  was  then  added  to  an  interpolated  surface  velocity  profile,  and 
this  velocity  profile  was  integrated  to  produce  a  height  profile  for  the  time  of  the  satellite 
pass  based  on  shipboard  data.  The  result  is  shown  in  Figure  6.3.  Even  when  an 
estimate  of  the  cyclostrophic  term  is  added,  the  height  difference  across  the  Gulf  Stream 
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is  about  12  cm  less  for  the  shipboard  determined  height  profile  than  it  is  for  the  height 
profile  measured  by  the  altimeter.  Interpolating  the  legs  4  and  7  velocity  profiles  tends 
to  smooth  out  features,  so  the  Gulf  Stream  and  any  other  features  will  tend  to  be  broader 
and  with  lower  peak  velocities  than  either  of  the  velocity  profiles  individually.  This 
would  explain  why  the  slope  of  the  height  curve  is  steeper  for  the  satellite  data  than  for 
the  shipboard  results,  and  why  the  curve  is  smoother  than  the  altimeter  curve,  but  does 
not  account  for  the  overall  difference  in  height  across  the  Gulf  Stream.  The  total  height 
difference  across  the  Gulf  Stream  is  should  be  the  same  for  both  curves  within  the 
expected  errors  for  the  two  techniques. 


1.4 

1.2 

1 


1       0.8 


Dynamic  Heights  at  time  of  day  107  GEOSAT  pass 


'I 
o 


0.6 


0.4 


0.2 


Dotted:  Geosat 

Solid:    Ship  geostrophic 

Dashed:   Ship  with  cyclostrophic  correction 


-0.2 1 


41  40  39  38  37  36  35  34  33  32 

Latitude 

Figure  63.  Sea  surface  height  for  the  day  107  GEOSAT  profile  compared  to  height  interpolated  between  ship's 
tracks  for  legs  4  and  7. 


One  source  of  error  is  the  time  difference  between  the  altimeter  measurements  and 
the  shipboard  measurements.  Interpolating  between  legs  4  and  7  to  the  time  of  the 
satellite  pass  is  an  attempt  to  reduce  this  problem.  The  height  difference  across  the  Gulf 
Stream  is  about  10  cm  larger  for  leg  4  than  for  leg  7.   Certainly  the  error  introduced  by 
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the  interpolation  is  less  than  ±  5  cm.  We  have  already  seen  how  the  ADCP  velocities  are 
prone  to  systematic  bias  errors.  Assuming  a  maximum  bias  of  no  more  than  2  cm/s  after 
balancing  the  leg  mass  transports  gives  a  maximum  error  of  ±  2.5  cm  in  dynamic  height 
across  the  Gulf  Stream. 

The  dynamic  height  computed  from  the  altimeter  data  is  also  subject  to  several 

sources  of  error: 

The  Kelly-Gille  model  for  the  mean  Gulf  Stream  is  based  on  processing  32  cycles 
of  the  satellite  between  November  1986  and  April  1988.  The  rms  height  error 
across  the  Gulf  Stream  due  to  the  height  variance  in  the  samples  used  to  compute 
the  mean  is  6  cm. 

The  particular  satellite  pass  is  subject  to  errors  caused  by  (1)  uncertainty  in  the 
measurement  of  the  travel  time  of  the  radar  signal,  and  (2)  uncertainty  in  the 
actual  orbit  position  of  the  satellite.  Height  measurement  errors  are  caused  by 
variability  of  the  moisture  content  of  the  troposphere  and  by  variations  in  the  free 
electron  content  in  the  ionosphere,  which  change  the  amount  of  refraction  of  the 
radar  signal,  and  by  errors  induced  by  variations  in  wave  height.  Orbit  errors  are 
due  to  variations  in  the  satellite  height  from  its  assumed  orbit  path.  The  expected 
error  in  height  difference  across  the  Gulf  Stream  due  to  these  effects  is  ±8  cm 
(Kelly,  et  al.  1990). 

The  total  rms  error  in  Ah  for  the  altimeter  curve  for  all  the  effects  described  above  is  ±10 

cm. 

The  height  profiles,  despite  the  difference  in  Ah  across  the  Gulf  Stream  are 
remarkably  similar,  with  several  smaller  scale  features  in  the  altimeter  profile  consistently 
reproduced  by  the  ship's  data,  despite  the  temporal  and  spatial  aliasing  problems.  More 
data  are  needed  from  future  cruises  to  determine  whether  better  agreement  is  achieved 
when  there  is  not  such  a  large  cyclostrophic  term  in  the  momentum  balance  and  when  the 
shipboard  and  satellite  data  are  taken  at  the  same  time. 
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Chapter  7 


Conclusions 


Velocities  measured  using  the  acoustic  doppler  current  profiler  have  been 
combined  with  geostrophic  velocity  profiles  generated  from  XBT  measurements  to 
estimate  the  flow  field  on  several  tracks  between  Bermuda  and  the  eastern  coast  of  the 
American  continent.  The  data  were  taken  during  two  separate  cruises,  the  first  in  April, 
1989,  and  the  second  in  December,  1989.  A  track  between  Bermuda  and  Cape  Cod  was 
chosen  to  correspond  to  an  ascending  subtrack  of  the  GEOSAT  altimeter,  allowing 
comparison  of  velocities  determined  from  shipboard  measurements  to  velocities 
determined  from  altimeter  data.   This  study  has  four  main  objectives: 

Demonstration  of  the  use  of  XBT's  and  the  ADCP  system  to  measure  upper  ocean 
transport  and  variability  over  large  areas  in  a  short  period  of  time. 

Use  of  repeated  measurements  along  the  Cape  Cod-Bermuda  track  in  order  to 
examine  short  term  variability  along  this  track  on  time  scales  of  order  one  to  two 
weeks. 
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Use  of  the  ADCP  velocities  to  compute  and  analyze  ageostrophic  terms  in  the 
horizontal  momentum  equations. 

Evaluation  of  the  use  of  satellite  altimetry  for  computing  sea  surface  velocities  by 
comparing  altimeter  results  to  shipboard  measurements  along  the  satellite  track. 

The  time  saving  advantage  of  the  ADCP-XBT  method  is  demonstrated  by  the  large 
areas  covered  during  the  two  cruises.  In  the  April  cruise,  for  example,  the  ship  covered 
a  distance  of  over  5000  kilometers  in  a  period  of  less  than  three  weeks  while  making 
continuous  ADCP  and  hourly  XBT  measurements.  This  allowed  measurements  of  the 
velocity  field  over  a  large  area  of  the  North  Atlantic  in  a  time  period  much  shorter  than 
that  required  for  hydrographic  measurements  with  CTD  stations.  The  disadvantages  of 
this  method  include  limitation  of  measurements  to  the  upper  ocean  and  errors  in  the 
geostrophic  velocity  profiles  introduced  by  the  inability  to  directly  measure  salinity. 

The  repeat  tracks  between  Bermuda  and  Cape  Cod  allowed  measurement  of  rapid 
variability  along  this  path.  Sources  of  variability  along  this  track  included  internal 
variability  in  the  Gulf  Stream,  which  caused  changes  in  transport  of  up  to  3  Sverdrups  in 
the  upper  765  meters  in  a  one  week  period,  and  a  rapid,  large  scale  barotropic  mode, 
which  caused  transport  changes  of  up  to  8  Sverdrups  in  one  week.  During  the  December 
cruise,  the  movement  of  a  cold  core  ring  onto  the  track  caused  a  change  in  transport  of 
8  Sverdrups  due  to  edge  effects  at  the  end  of  the  track. 

Measurements  of  horizontal  velocity  components  along  the  repeated  track  using 
the  ADCP  system  allowed  computation  of  both  temporal  and  spatial  derivatives  along  the 
track.  This  allowed  us  to  measure  the  size  of  some  of  the  ageostrophic  terms  and  to 
estimate  the  size  of  the  errors  introduced  in  using  the  geostrophic  approximation. 
Measurement  of  the  non-geostrophic  terms  indicates  that  the  geostrophic  approximation 
is  reasonably  good  for  computing  the  along-stream  velocity  (errors  <,  10%)  but  that 
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geostrophy  is  a  poor  approximation  for  deterniining  the  cross-stream  velocity  component 
(errors  20-35%). 

Sea  surface  velocity  measurements  made  from  the  GEOSAT  altimeter,  when 
compared  to  measurements  made  from  the  ship  along  the  same  track,  show  a  consistent 
picture  of  the  size,  location,  and  motion  of  the  Gulf  Stream.  Smaller  scale  features  in  the 
flow  field  are  also  reproduced  consistently,  despite  some  temporal  and  spatial  aliasing. 
These  measurements  were  taken  in  the  trough  of  a  large  Gulf  Stream  meander,  which 
introduced  a  significant  cyclostrophic  term  into  the  momentum  balance.  This  term  was 
estimated  using  satellite  imagery  to  compute  the  radius  of  curvature  of  the  Gulf  Stream, 
and  subtracted  from  the  shipboard  data.  The  geostrophic  GEOSAT  velocities  in  the  Gulf 
Stream  exceeded  the  shipboard  measurements  by  about  10%,  even  after  the  cyclostrophic 
term  was  removed.  The  difference  between  the  shipboard  and  altimeter  estimates  of 
dynamic  height  are  attributed  to  errors  in  the  altimeter  height  profile,  caused  by 
limitations  of  the  model  used  to  remove  the  mean  Gulf  Stream  from  the  geoid,  and  by 
inaccuracies  in  computing  the  radar  signal  path  and  the  orbit  path  of  the  satellite. 
Additional  comparisons  of  shipboard  and  altimeter  measurements  for  different  conditions 
in  the  Gulf  Stream  are  needed  to  further  evaluate  the  accuracy  of  satellite  altimetry  for 
determining  sea  surface  current  velocities. 
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